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This  contract  was  initiated  to  determine  the  aeroelastic  stability 
limits  of  articulated  and  unarticulated  helicopter  rotor  systems  at 
high  forward  speeds.  The  four  primary  modes  of  aeroelastic  instability 
(classical  flutter,  stall  flutter,  torsional  divergence,  and  flapping 
or  flatwise  bending  instability)  were  investigated.  The  possibility 
of  a  flap-lag  instability  suggested  by  Dr.  Maurice  X.  Young  of  the 
Vertol  Division,  The  Boeing  Company,  was  inves  igated  as  a  special  case 
of  flapping  instability. 


The  results  are  published  as  a  five-volume  oet;  the  subject  of  each 
volume  is  as  follows: 


Volume  I 
Volume  II 
Volume  III 
Volume  IV 
Volume  V 


Equations  of  Motion 
Classical  Flutter 
Stall  Flutter 
Torsional  Divergence 
Flapping  Instability 


These  reports  have  been  reviewed  by  the  U.  S.  Army  Aviation  Materiel 
Laboratories.  These  reports,  which  are  published  for  the  exchange  of 
information  and  the  stimulation  of  ideas,  are  considered  to  be  tech¬ 
nically  sound  with  regard  to  technical  approach,  results,  conclusions, 
and  amended  parameter  ranges  for  accurate  usage. 
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SUMMARY 


The  purposes  of  this  research  program  were  to  extend  or  develop  analytical 
methods  for  determining  rotor  Linde  aeroelastie  stability  limits  and  to 
perform  stability  calculations  over  a  range  of  design  and  operating  vari¬ 
ables  for  articulated  and  nonarticulated  configurations.  The  usefulness 
of  simpler  analytical  methods  was  investigated  by  Comparing  results  with 
operating  boundaries  from  the  more  elaborate  analysis. 

The  analytical  study  in  this  volume  was  carried  out  to  determine  the 
susceptibility  of  helicopter  rotor  blades  to  a  stall  flutter  instability. 
This  analysis  was  based  on  the  use  of  unsteady  aerodynamic  data  previously 
obtained  by  Sikorsky  Aircraft,  for  an  HACA  0012  airfoil  oscillating  in  pitch 
about  its  quarter-chord  over  a  wide  range  of  values  of  incidence  angle, 
oscillatory  frequency,  amplitude  of  motion,  and  free-stream  velocity. 

These  data  were  originally  available  in  the  form  of  moment  coefficient- 
incidence  angle  loops,  and  a  twofold  task  was  performed  i.i  carrying  out 
the  present  study.  First,  it  was  necessary  to  convert  the  moment,  coeffi¬ 
cient  data  to  an  aerodynamic  damping  parameter  form.  This  was  accomplished 
by  integrating  the  moment  over  one  cycle  of  motion  to  yield  the  aerodynamic 
work  per  cycle,  and  this  in  turn  was  multiplied  by  appropriate  conversion 
factors  to  produce  the  desired  two-dimensional  aerodynamic  damping. 

Second,  it  was  necessary  to  apply  these  two-dimensional  results  to  a  heli¬ 
copter  rotor  to  evaluate  the  weighted  three-dimensional  damping  at  each 
azimuth  station,  and  to  interpret  the  implications  of  any  predicted  region 
of  instability. 

It  was  found  that  under  certain  combinations  of  forward  speed  and  high 
disc  loading,  a  helicopter  rotor  blade  could  encounter  regions  of  negative 
torsional  aerodynamic  damping  while  operating  as  a  retreating  blade.  The 
extent  of  the  negative  damping  regime  was  sensitive  to  change  in  'ncidenee 
angle  distribution,  and  for  some  parameter  combinations  it  could  produce 
two  or  three  cycles  of  unstable  torsional  motion  per  revolution, 

Generally  good  agreement  was  found  to  exist  between  the  results  of  the 
present  study  and  those  of  other  investigators  at  ail  stages  in  the  devel¬ 
opment;  e.g.  ,  the  two-dimensional  damping  data  were  in  good  agreement  with 
the  results  of  Reference  1  and  the  final  three-dimensional  damping  vari¬ 
ation  with  azimuth  over  the  rotor  disc  was  in  qualitatively  good  agreement 
with  the  recent  work  of  Reference  2. 

The  stall  flutter  analysis  was  used  in  conjunction  with  the  blade  motion 
solution  of  Volume  I  to  provide  flight  condition  boundaries  for  stall 
flutter  intensity. 
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The  investigation  presented  in  this  volume  is  part  cl'  an  extensive  study, 
which  3 s  presented  in  five  'volumes-  The  work  vbo  performed  under  Contract 
PA  ^-17T-AMC-33C(T)  witli  the  l) .  G .  Army  Aviation  Materiel  Laboratories, 
Fort  Fustic,  Virginia.  The  program  was  monitored  1'or  U5AAVLA13S  by 
Mr.  Joseph  Me Garvey . 

The  rotor  blade  stall  Flutter  analysis  presented  in  this  volume  is  the 
work  of  Mr.  Franklin  0.  Carta  of  the  United  Aircraft  Research  Laboratories, 
The  test  data  used  in  the  analysis  are  from  a  wind  tunnel  test  which  was 
performed  previously  under  Sikorsky  Aircraft  sponsorship. 

The  subheading  entitled  "Stall  Flutter  Flight  Condition  Boundaries"  is 
the  result  of  work  done  at  Sikorsky  Aircraft  by  Mr.  Charles  F-  Niebanck. 

The  information  under  that  subheading  is  the  result  of  the  application  of 
the  stall  flutter  analysis  of  this  volume  in  conjunction  with  the  Extended 
Normal  Mode  Transient  .Analysis  of  Volume  1. 

Volume  1  of  this  report  contains  the  development,  of  the  differential 
equations  of  motion  of  an  elastic  rotor  blade  with  chordwise  mass  un¬ 
balance,  . 

Volume  11  presents  a  linearized  discrete  azimuth  classical  flutter  analysis 
for  rotor  blades,  witli  an  appropriate  parameter  variation  study,  a  com¬ 
parison  with  test  data,  and  a  comparison  with  results  calculated  by  using 
the  method  of  Volume  I. 

Volume  IV  contains  the  results  of  a  study  of  static  torsional  divergence. 

A  set  of  design  charts  and  the  effects  of  a  range  of  parameter  variations 
are  presented.  The  results  of  the  stati  •  divergence  calculation  are  com¬ 
pared  with  results  calculated  by  using  the  method  of  Volume  I. 

Volume  V  presents  the  results  of  a  study  of  flapping  and  coupled  flap-lag 
instabilities.  The  results  cf  a  parametric  study  based  on  a  single  degree- 
of-freedom  flapping  or  flatwise  bending  analysis  are  presented.  Compari¬ 
sons  are  made  with  results  from  the  more  elaborate  method  ol  Volume  I. 

The  results  obtained  by  using  the  method  of  Volume  I  to  determine  the 
coupled  flap-lag  response  of  a  rotor  to  a  number  of  sudden  control  changes 
are  presented.. 
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INTRODUCTION 


Recently,  as  a  result  of  tlie  combined  requirements  of  higher  forward  speeds 
and  greater  blade  loadings  for  helicopters ,  there  has  been  a  renewed  in¬ 
terest  in  the  study  of  stall-induced  oscillatory  blade  instabilities,  with 
particular  emphasis  on  the  effects  of  these  extreme  requirements  on  the 
stall  flutter  phenomenon.  This  is  evidenced  by  the  work  reported  in 
References  2,  3,  and  b  ,  to  name  only  a  few.  However,  with  the  exception 
of  Reference  2,  all  of  these  studies  had  one  major  deficiency  in  common  — 
a  dearth  of  valid,  unsteady  aerodynamic  data  encompassing  a  substantial 
portion  of  the  stall  flow  regime.  3n  fact,  even  in  Reference  2  the  data 
were  synthesized  and  idealized  from  a  number  of  sources. 

These  increased  performance  requirements  were  foreseen  some  time  ago,  and 
Sikorsky  Aircraft  obtained  a  large  quantity  of  two-dimensional  unsteady 
aerodynamic  test  data  from  an  isolated  NAOA  0012  airfoil  section,  which 
was  oscillated  in  pitch  about  its  25%  chord.  These  data  are  utilized  in 
the  stall  flutter  stability  analysis  described  in  this  volume. 

Similar  data  were  obtained  from  Reference  1.  These  data  were  converted  to 
a  form  comparable  with  the  Sikorsky  data  and  were  used  in  a  similar  sta¬ 
bility  analysis.  The  two  sets  of  results  showed  good  qualitative  agreement. 

The  stall  flutter  stability  analysis  was  used  to  provide  stall  flutter 
boundaries  for  a  typical  rotor.  The  boundaries  indicate  flight  conditions 
for  which  stall  flutter  may  be  particularly  troublesome. 


DESCRIPTION  OF  TEST  EQUIPMENT,  MODEL  INSTALLATION ,  AND  PROCEDURE 


TEST  EQUIPMENT1  AND  MODEL 


The  experimental  program  was  conducted  in  the  two-dimensional  channel  of 
the  UAC  8-foot  main  wind  tunnel,  which  has  an  atmospher ic ,  single-return, 
closed-throat  circuit  with  an  octagonal  throat.  Test  section  Mach  number 
is  variable  below  M  =  1.0,  and  the  Reynolds  number  for  the  2-foot  chord 
model  used  in  the  test  program  was  approximately  Rc  -  1.4  M  x  10?.  Figure 
1  shows  schematically  the  mechanical  system  used  in  the  dynamic  tests  to 
obtain  unsteady  aerodynamic  data  for  the  oscillating  airfoil. 

The  NACA  0012  model  tested  during  this  program  consisted  of  a  balsa  wood 
core  supported  by  chordwise  aluminum  stiffeners  and  wrapped  in  three  layers 
of  fiber  glass.  The  model,  which  had  a  span  of  33  inches  and  a  chord  of 
24.5  inches,  is  shown  in  Figure  2.  Thirteen  pairs  of  differential  pressure 
orifices  were  located  on  the  top  and  bottom  of  the  model  along  a  chordwise 
line  which  was  11.5  inches  inboard  from  the  tip.  The  orifices  were  connec¬ 
ted  with  internal  tubing  to  13  miniature  differential  pressure  pickups  of 
the  variable  reluctance  type  which  were  enclosed  in  a  cavity  in  the  model 
tip.  In  this  manner,  pressure  acting  on  the  airfoil  surface  at  a  given 
chordwise  station  was  converted  to  an  electrical  signal  which  was  trans¬ 
mitted  to  a  recording  oscillograph.  Representative  oscillograph  traces  of 
the  pressure  transducer  signals  are  presented  in  Figure  3.  In  addition,  a 
linear  transformer  was  attached  to  the  end  ol'  the  shaft,  and  its  output  was 
used  to  record  variations  in  angle  of  attack. 

TEST  PROCEDURE 

In  general,  the  following  test  procedure  was  employed.  A  set  of  eccentric 
cams  was  installed  which  fixed  the  amplitude  of  motion,  a  .  After  a  pre¬ 
liminary  shake  test,  all  recording  instruments  were  zeroed  and  the  tunnel 
was  started  and  brought  up  to  its  operating  Mach  number  M  .  At  the  pre¬ 
scribed  Mach  number,  the  mean  angle  of  attack,  ,  was  varied  over  a 

range  of  values,  and  for  each  value  of  the  frequency,  f  ,  was  varied 

over  a  range  of  values.  After  data  had  been  recorded  at  all  of  the  desired 
combinations  of  and  f  ,  the  Mach  number  was  changed  and  records  were 

taken  for  a  new  set  of  values  of  av  and  <  .  Finally,  after  all 
desired  values  of  M  had  been  obtained,  the  cams  were  changed  to  obtain 
a  new  value  of  a  and  the  entire  process  was  repeated. 

The  oscillograph  records  of  the  steady  and  unsteady  pressure  response  were 
taken  for  most  combinations  of  the  following  parameter  values: 

Amplitude  of  motion,  Q  -  +  4,  6,  8  deg 

Mach  number ,  M  0.2,  0.3,  0*4 

Mean  angle  of  attack,  =  0  to  33  deg  in  3-deg  increments 

Frequency,  1  =  0.0,  0.5,  1,  2,  4,  8,  12,  16  cps 


In  the  interest  of  expediency  or  safety,  certain  combinations  were  omitted. 
For  example,  it  was  decided  to  limit  the  maximum  frequency  to  12  cps  for 
mean  angles  of  attack  greater  than  2k  degrees  to  ensure  the  structural 
integrity  of  the  system,  particularly  at  the  higher  Mach  numbers. 
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DATA  REDUCTION 


The  oscillograph  records  were  read,  and  the  measured  data  were  tabulated. 

A  digital  computer  was  then  used  to  reduce  these  data  to  obtain,  for  each 
run,  the  chordwise  differential  pressure  distribution  as  well  as  normal 
force  and  moment  as  functions  of  instantaneous  angle  of  attack.  The  mag¬ 
nitude  of  the  corrections  to  the  measured  pressures  for  the  influence  of 
tunnel  wall  constraint  under  steady -state  conditions  was  estimates  to  be 
less  than  3!t  of  the  full-scale  values  in  the  linear  range;  therefore,  no 
corrections  to  the  data  were  made.  The  unsteady  wall  interference  effects 
were  small  compared  to  the  gross  unsteady  parameters  being  measured  and 
were  also  neglected  (Reference  5). 

FRESSURE  DISTRIBUTIONS 


The  output  of  each  pressure  transducer  was  directly  proportional  to  the 
local,  instantaneous  pressure  difference. 


Ap*(  x,  t)  =  Ap*(x)  = 


upper 


n 

tower 


(1) 


where  the  asterisk  denotes  a  complex  quantity  and  the  superscript  bar^ 
denotes  an  amplitude  function.  The  pressure  difference  amplitude  Ap 
is  expressed  as  a  complex  number, 

Ap  *-  ApR  +  iApj  (f) 


to  account  for  possible  phase  differences  between  the  local  pressure 
response  and  the  motion.  Ac  shown  in  Figure  3,  the  peak-to-peak  difference 
in  the  oscillograph  trace  was  proportional  to  the  absolute  magnitude  of  the 
complex  pressure  difference  given  in  Eq.  (2), 
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Scxne  representative  results  of  these  measurements  a re  shown  in  Figures  ^ 
and  5,  in  which  the  absolute  magnitude  of  the  pressure  difference  from 
Eq.  (3)  has  been  divided  by  the  dynamic  pressure, 

q  -  y  pu2  (M 
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und  the  resulting  pressure  coefficient,  lAp/ql  ,  has  been  plotted  as  a 
function  of  dimensionless  chordal  .-.ration.  In  Figure  ii,  the  experimentally 
determined  values  of  this  dimensionless  amplitude  of  the  differential 
pressure  coefficient  ore  presented  for  a  range  of  mean  angles  of  attack 
(  Qm  =  0  to  32  degrees  in  3-degree  increments)  at  u  fixed  frequency  of  it 
cps.  In  Figure  5, a  similar-  set  of  results  is  presented  for  a  range  of 
frequencies  (f  =  it ,  8,  IF,  16  ops)  for  a  given  set  of  mean  angles  of  attack 
( aM  =0,  3,  6,  9  degrees).  In  both  figures,  the  solid  curve  represents  the 
pressure  coefficient  predicted  from  classical  unsteady  potential  flow 
theory.  In  the  present  case  the  theoretical  pressure  distribution  formula 
was  obtained  from  Reference  6,  which  is  based  on  the  fundamental  work  of 
Theodorsen  (Reference  ?)•  The  pitching  results  in  Reference  6  have  been 
derived  for  a  pivot  axis  at  the  2.'?%  chord;  for  convenience,  the  real  and 
imaginary  parts  are  reproduced  below,  with  the  notation  revised  slightly 
to  conform  to  present,  usage. 


In  these  formulas,  x  is  the  dimensionless  chordwise  position  relative  to 
the  raidchord; in  semichords,  and  k  is  the  reduced  frequency  parameter. 


and  F  and  G  are  the  real  and  imaginary  parts,  respectively,  of  the 
Theodorsen  circulation  function,  C(k)  -F(k)  +  iG(k). 

It  is  seen,  from  Figure  ,  that  at  low  to  moderate  incidence  angles,  the 
experimental  values  are  in  good  agreement  with  theory;  even  at  high  in¬ 
cidence  angles,  the  leading-edge  peak  is  still  maintained,  although  the 
agreement  between  theory  and  experiment  has  deteriorated.  Finally,  at 
extremely  high  incidence  angles,  even  the  leading-edge  experimental  values 
no  longer  agree  with  theory .  The  effects  of  increasing:  frequency  at  low 
incidence  angle  are  shown  5:i  Figure  5, und  it  is  seen  that  in  all  four  parts 
of  this  figure  the  agreement  between  theory  ana  experiment  is  excellent. 


NORMAL  FORCE  AND  MOMENT  AMPLITUDE; 


On  pages  39 ^  and  395  of  Reference  8  are  found  the  equations  for  normal 
force  and  moment  about  the  pivot  axis  in  terms  of  the  pressure  distribu¬ 
tions  , 


hi*  = 
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or,  in  coefficient  form, 
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Pressure  data  obtained  during  the  experimental  program  were  reduced  to 
coefficient  form  by  use  of  Eqs  .  (10)  and  (ll).  The  results  of  these  in¬ 
tegrations  form  the  basis  for  the  stability  analysis  to  follow.  However, 
before  proceeding  with  this  analysis,  consideration  will  be  given  to  Figure 
6  in  which  the  amplitudes  of  Cn  and  Cm  have  been  plotted  versus 
reduced  frequency,  k  ,  over  a  small  range  of  mean  angles  of  attack. 
Included  in  both  of  these  figures  are  the  theoretical  variations  of  these 
functions  as  predicted  by  the  potential  flow  analysis  (Reference  7).  The 
analytic  cl  expressions  used  here  may  be  obtained  from  Reference  8,  and  the 
explicit  relationships  for  unsteady  moment  cueff  ±C  I  cliti  will  be  considered 
in  great  detail  in  a  subsequent,  section  of  this  volume.  For  the  moment, 
though,  it  is  sufficient  to  note  that  the  general  trends  of  the  experi¬ 
mental  data  are  in  good  agreement  with  the  theory  for  mean  incidence  angles 
up  to  Qm  =  9  degrees;  in  fact,  the  data  for  zero  mean  incidence  angle  are 
in  excellent  agreement  witn  theory. 

NORMAL  FORCE  AND  MOMENT  HYSTERESIS  LOOPS 

When  the  instantaneous  normal  force  or  moment  coefficient  is  plotted  versus 
angle  of  attack,  the  resulting  exosed  curve  surrounds  an  area  which,  in  the 
case  of  the  moment,  is  representative  of  the  energy  absorbed  or  dissipated. 
This  is  found  to  be  the  case  in  both  classical  and  separated  flows.  Con¬ 
sider  the  case  depicted  in  Figure  7-  The  two  upper  curves  represent  hypo¬ 
thetical  sinusoidal  variations  of  angle  of  attack  and  either  normal  force 
or  moment,  both  expressed  in  arbitrary  units,  as  a  function  of  time.  In 
the  example  shown,  the  normal  force  or  moment  leads  the  angle  of  attack  by 
one-eighth  of  a  cycle  or  15  degrees  in  phase.  In  the  bottom  portion  of 
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the  figure ,  the  time  variable  has  been  eliminated  and  the  normal  force  or 
moment  has  been  plotted  directly  as  a  function  of  angle  of  attack.  The 
arrows  denote  the  direction  of  increasing  time.  This  figure  shows  that  a 
phase  shift  between  force  or  moment  and  motion  produces  a  ] oop  which  en¬ 
closes  a  finite  area.  The  same  effect  will  be  found  to  exist  in  the  case 
of  a  nonsi nusoidal  force-motion-time  relationship  caused  by  the  presence 
of  separated  flew.  in  classical  or  potential  flow,  the  closed  contour  will 
be  elliptic.il;  whereas  in  separated  flow,  the  contour  will  be  distorted. 

Normal  Force  Loops 

A  few  representative  unsteady  normal  force  coefficient  loops  are  presented 
in  Figure  8  for  a  constant  amplitude  of  <5  ~  6  degrees  and  a  constant 

Mach  number  of  M  =  0.3.  In  this  figure  the  experimentally  determined  un¬ 
steady  normal  force  coefficients  are  plotted  versus  instantaneous  incidence 
angle,  with  the  superimposed  arrowheads  indicating  the  direction  of  in¬ 
creasing  time.  The  solid  lines  represent  the  unsteady  data, and  the  dashed 
lines  represent  the  steady- state  characteristics  (also  obtained  from 
pressure  readings  in  this  test).  The  arrows  indicate  the  direction  of  in¬ 
creasing  time.  Three  of  the  inset  figures  were  for  a  constant  frequency 
of  4  cps  and  serve  to  snow  the  effects  of  varying  mean  incidence  angle, 
from  qm  =  6  to  12  to  18  degrees.  The  two  right-hand  inset  figures  are 
both  for  a  mean  incidence  angle  of  <3m  =  12  degrees  and  illustrate  the 
effect  of  a  change  in  frequency  from  4  cps  to  16  cps .  It  is  clear  that  the 
increase  in  mean  incidence  angle  to  values  greater  than  the  steady-state 
stall  angle  has  a  rather  profound  effect  on  the  dynamic  force  response  of 
the  oscillating  airfoil.  It  is  also  clear  that  an  increase  in  frequency 
produces  a  radical  change  in  the  dynamic  stalling  behavior  of  the  airfoil. 
Specifically,  at  low  frequency  the  dynamic  force  response  reaches  its  peak 
value  just  before  the  maximum  incidence  angle  is  reached;  it  then  drops 
precipitously  to  a  value  far  below  the  steady-state  stall  value  and  remains 
tntre  for  almost  the  entire  region  of  decreasing  incidence.  In  contrast  to 
this  behavior,  the  effect  of  nigh  frequency  is  to  maintain  a  nearly  ellip¬ 
tical  response  loop,  even  for  incidence  angles  beyond  stall,  over  the 
entire  range  of  instantaneous  incidence  angle.  The  further  significance 
of  this  behavior  will  be  explored  at  length  below  and  in  subsequent  sections 
in  connection  with  the  moment  coefficient  loops  and  their  stability  impli¬ 
cations  . 

Moment  Loops 

Some  representative  unsteady  moment  coefficient  loops  are  presented  in 
Figures  9  through  11.  These  loops  are  the  basis  for  the  stability  analysis 
which  follows;  hence, they  are  shown  in  greater  detail  than  the  normal  force 
jooris.  In  these  figures,  an  effort  has  been  made  to  illustrate  the  vari¬ 
ations  in  the  moment  hysteresis  loops  for  each  of  the  variable  test  param¬ 
eters.  Thus,  the  effect  of  varying  mean  incidence  angle,  Cif/  ,  is  shewn 
in  Figure  9;  the  effect  of  varying  frequency,  f  ,  is  shown  in  Figure  10; 
and  the  effect  of  varying  Mach  number,  M  ,  and  torsional  amplitude,  q  > 
is  shown  in  Figure  11.  In  all  of  these  figures, the  superimposed  arrowheads 
indicate  ■  direction  of  increasing  tiiae. 
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Figure  9  illustrates  the  variation  in  moment  coefficient  with  incidence 
angle  at  the  fixed  conditions  M  =  0.2,  a  =6  degrees  and  f  =  U  cps . 

The  solid  curves  are  associated  with  the  unsteady  motion  presently  under 
consideration,  and  the  dashed  curve  represents  the  steady-state  moment 
(from  balance  data)  for  the  same  Mach  number  and  Reynolds  number.  (Once 
again, the  arrows  indicate  the  direction  of  increasing  time.)  It  can  be 
seen  that  all  three  moment  loops  generally  follow  the  steady-state  curve, 
and  it  will  be  shown  presently  in  Figure  10  that  the  loop  for  <2^  =  0 
degrees  is  in  good  agreement  with  the  results  obtained  from  potential  flow 
theory.  (This  is  also  evidenced  by  the  proximity  of  the  data  point  for 
Ciu  -  0  degrees  to  the  theoretical  curve  in  Figure  6  for  a  reduced  fre¬ 
quency  value  of  k  =  0.112.)  The  character  oi  tne  hysteresis  loop  changes 
radically  as  aM  increases,  and  the  moment  loop  for  aM  =  12  degrt  s  dis¬ 
plays  the  characteristic  crossover  behavior  of  the  unsteady  moment  in  the 
neighborhood  of  the  stalling  angle,  Tnis  separation- induced  distortion  of 
the  moment  hysteresis  loop  has  been  the  subject  of  a  number  of  early  in¬ 
vestigations,  including  those  reported  in  References  1  and  9  and,  hence, 
will  not  be  discussed  at  length  herein.  However,  it  will  be  shown  later 
in  this  report  that  the  system  stability  is  strongly  dependent  on  the 
direction  in  which  the  area  of  the  moment  loop  is  enclosed;  in  particular, 
it.  will  be  shown  that  a  counter  clockwise  enclosure  is  stable,  whereas  a 
clockwise  enclosure  is  unstable.  Thus,  it  is  seen  that  in  some  mean  in¬ 
cidence  range  containing  =  i2  degrees,  there  may  be  some  potentially 
unstable  regions  which  could  affect  the  stability  of  rotor  systems. 

Finally,  for  a  mean  incidence  angle  of  2k  degrees,  the  moment  loop  once 
again  implies  a  stable  motion. 

Figure  10  shows  the  effects  of  increasing  - iie  frequency  from  U  cps  to  16 
cps  for  two  mean  angles  of  incidence.  Also  in  Figure  10,  the  experimental 
results  are  shown  as  solid  curves,  and  the  results  from  unsteady  potential 
flow  theory  are  shown  as  dashed  curves.  (A  brief  discussion  of  these 
theoretical  predictions  is  presented  later  in  this  volume  in  Appendix  I). 

It  is  seen  that  except  for  a  slight  upward  displacement  of  the  experimental 
results  relative  to  the  theoretical  curves,  the  two  are  in  excellent  agree¬ 
ment.  The  right  side  of  this  figure  shows  that  for  CiM  =  15  degrees,  the 
effect  of  an  increase  in  frequency  (and,  hence,  an  increase  in  reduced 
frequency)  appears  to  be  stabilizing,  even  though  the  incidence  angle  is 
considerably  greater  than  the  stalling  angle  over  much  of  the  range.  This 
is  in  accord  with  the  implied  results  in  References  10  and  IT,  as  pointed 
out  and  further  amplified  by  Reference  1,  and  also  as  reported  in  Reference 
12.  However,  in  the  left  column  of  Figure  11,  it  appears  that  an  increase 
in  reduced  frequency  (caused  by  a  decrease  in  Mach  number  at  constant  fre¬ 
quency,  f  =  k  cps)  yields  a  contradictory  result  in  that  the  stable 
closed  loops  for  k  =  0.056  and  0.075  are  replaced  by  a  marginally  stable 
crossed  loop  for  k  =  0.112.  It  is  believed  that  this  behavior  is  char¬ 

acteristic  of  very  small  values  of  k  ,  and  it  will  be  seen  later  in  this 
volume  that  all  of  these  results  are  self-consistent.  Finally,  the  right 
column  of  Figure  11  shows  the  effect  of  increasing  the  torsional  amplitude. 
Superimposed  on  each  curve  in  Figure  II  is  the  steady-state  moment  vari¬ 
ation,  represented  by  the  dashed  line.  Once  again,  the  results  appear  to 
be  self-consistent. 


This  concludes  the  review  of  the  original  test  program.  In  the  next  sec¬ 
tion,  the  analytical  expressions  necessary  fei  the  stability  investigation 
will  be  developed;  following  this,  the  application  of  this  analysis  to  the 
data  will  be  discussed.  The  last  section  of  the  report  will  deal  with  the 
stability  characteristics  of  a  typical  rotor  under  various  loading  con¬ 
ditions  . 


ANALYSIS 


TWO-DIMENSIONAL  WORK  PER  CYCLE 


T'ae  differential  work  done  by  the  aerodynamic  moment  during,  the  course  of 
the  torsional  motion  is  obtained  by  computing  the  product  of  the  in-phase 
components  of  moment  and  differential  twist,  or 

tfw  =  M  nd  aH  (12) 


where  Mr  and  Qr  are  the  real  parts  of  these  quantities.  Hence,  the 
work  per  cycle  of  motion  is  obtained  by  integrating  Eq.  (12)  over  one 
cycle,  or 


W 


(13) 


This  may  be  rendered  dimensionless  by  dividing  both  sides  of  the  equation 
by  ~p U2(2b)Z  >  which  yields  the  equation  for  the  work  coefficient  in 
terms  of  the  integral  of  the  moment  coefficient,  as  follows: 


Cw 
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To  evaluate  the  integral  in  Eq.  (lU),  it  is  necessary  to  introduce  the 
quantities  Cm»  and  ,  The  derivation  of  these  functions  is  presented 

in  Appendix  I,  which  lias  been  included  in  the  present  volume  for  complete¬ 
ness.  From  Eq.  (71*)  ,  after  expanding  the  exponential  function  in  sines 
and  cosines,  the  real  part  of  the  moment  coefficient  may  be  obtained  in 
the  form 


^ij  "  Cl.  +  C,.  cos  u>1  -  Cjy,  sin  out 
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The  differential  of  Eq.  (80)  is 


and  after  Eqs .  (15)  and  (lo)  are  substituted  .into  Eq.  (ll)»  the  result  is 
Cw  5  -  f2ir  j  cm  +  C  cos  tuf  -  C  sin  cut  a  sm  cut  d  M )  (17) 
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(The  integration  range,  0<  <u<  27T  ,  is  equivalent-  to  one  complete  cycle 

of  motion.)  After  the  integrations  indicated  in  Eq .  (17)  are  performed, 
it  is  found  that  the  term  involving  the  mean  moment  vanishes  as  well  as 
the  term  containing  the  real  part  of  the  unsteady  moment;  the  final  result 
for  the  theoretical  work  coefficient  is  given  by 

C--nSCu  (18) 
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This  is  the  work  done  by  the  air  on  the  airfoil;  hence,  a  positive  value 
of  Cw  will  indicate  an  unstable  motion,  since  this  implies  a  net  energy 
exchange  from  the  surrounding  medium  to  the  airfoil,  whereas  a  negative 
value  of  Cw  will  indicate  a  stable,  or  dumped  motion. 

TWO -AMD  THREE-DIMENS JOKAL  AERODYNAMIC  DAMPING 

Before  proceeding  with  the  analysis  of  the  aerodynamic  damping  of  the  sys¬ 
tem,  it  is  useful  to  review  briefly  the  behavior  of  a  linear,  damped,  tor¬ 
sional  system  such  as  the  one  described  by  the  differential  equation 

I  a*  ♦  c  a*  +  xa’  -  0  (19) 

where  I  is  the  inertia,  C  is  the  damping,  and  k  is  the  stiffness  of 
the  system.  If  the  motion  is  essentially  sinusoidal  (i.e.,  only  slightly 
damped  and,  hence  very  nearly  a  constant  amplitude  sinusoid)  then  Eq.  (69) 
of  Appendix  I  is  a  solution.  The  equation  becomes 

(-<uzI  +  i<uc  +  <  )  a  -  0  (20) 


It  is  seen  that  the  damping  coefficient  is  contained  in  the  imaginary  part 
of  this  expression,  and  it  may  be  assumed  that  the  equivalent  damping  terms 
for  any  similar  linear  system  will  also  be  contained  in  the  imaginary  part 
of  the  differential  equation  solution. 

Eq.  (19)  represents  a  syst  em  oscillating  in  torsion  in  a  vacuum.  If  the 
same  system  were  to  oscillate  in  torsion  in  a  moving  airstream,  the  right- 
hand  side  of  the  equation  would  no  longer  be  zero  but  would  represent  the 
unsteady  moment  imposed  by  the  moving  air  on  the  body,  or 

la  +  Ca  +  xa  '-My  (2l) 
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It  is  shown  in  References  8  and  13  that  in  the  case  of  a  single-degree-of- 
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torsional  displacement  and  its 
pressed  in  the  general  form 


"the  unsteady  inoiuent  is  a  function  of  the 
first  two  time  derivatives  and  may  be  ex- 
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In  the  results  cited  in  References  8  and  13  for  an  isolated  airfoil  oscil¬ 
lating  in  an  incompressible,  potential  flow  (based  on  Theodorsen's  theory, 
Reference  7),  the  component  M>  is  a  pure  real  quantity,  whereas  both 
and  are  complex;  hence. 


M„*  =  M 


2R  4 


i  M, 


21 


<  -  MJR  ♦ 


iM 


31 


(23) 


After  Eqs .  (22)  and  (23)  are  substituted  into  Eq.  (2l)  and  the  terms  are 
rearranged  slightly,  the  result  is 

(I-M,)d%(C-M2R-  iM  2iy*  (k-Mw-  IM  3I)a  =0  (  2h ) 


Once  again,  if  the  damping  is  sufficiently  small,  the  motion  will  be  nearly 
sinusoidal  and  Eq.  (69)  represents  a  solution,  whereupon  Eq.  (2^)  becomes 


[  “2  (I  M.) 4  i"  (c  -  M2R  iM3i) 


a  T  0 


(25) 


After  collecting  real  and  imaginary  parts, 

[{" 41,2  (*  ‘  Mi)  4  ujM2i+  k  -  i  (C  -  M2R)-  M3i|]  3=0  (  ro) 


As  in  the  case  of  the  system  oscillating  in  vacuum,  the  imaginary  part  of 
Eq.  (26)  will  represent  the  total  damping  of  the  system,  as  follows: 

total  damping  =  u>  (0  -  M2R)  '  ^31  (27) 

Since  the  quantity  cue  constitutes  the  system  damping  in  the  absence  of  a 
moving  airstream,  then  the  remainder  of  Eq.  (27)  must  be  the  effective 
aerodynamic  damping  of  the  system,  as  follows: 

£=-WM2r-M3i  (28) 
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The  symbol  (  denotes  the  (dimensional)  aerodynamic  damping  parameter  of 
a  system  executing  a  si ngle-degree-of-freedom  torsional  motion.  Return 
now  to  Eq .  (22).  It  will  be  assumed  that  the  left-hand  side  of  this  equa¬ 
tion  is  a  time  dependent,  sinusoidal  function  having  the  form 


(29) 


where 


After  Eqs.  (29),  (30),  and  (69)  are  substituted  into  Eq.  (22)  and  the  ex- 
poential  factor  elu,t  is  cancelled,  the  result  is 

Mun  +  *  Mm  -  -  ui2  M,  -  Mjr]  +  i  [cu  M2R  4  W3I  ]Ja  (  31 ) 


When  real  and  imaginary  parts  are  equated,  it  is  seen  that 


Mul  : 


(wM2r  +  M31. 


a 


(32) 


A  comparison  of  Eqs.  (28)  and  (32)  shows  that  the  aerodynamic  damping 
parameter,  £  ,  is  equal  to  the  negative  of  the  derivative  of  the  imaginary 
component  of  the  unsteady  moment  with  respect  to  the  amplitude  of  motion. 


d  Mui 

da 


-  [wM?r+  M3I] 


(33) 


It  is  convenient  at  this  point  to  rewrite  Eq.  (33)  in  dimensionless  coef¬ 
ficient  form,  as  previously  done  in  Eq .  (1^).  Hence,  after  defining  the 
dimensionless  two-dimensional  aerodynamic  damping  parameter  in  pitch  by 
the  symbol 


t?  r  _ i _ 

~q2  (l/2)^uZ(2b)2 


(3V 


then  Eq.  (33)  may  be  rewritten  as 


q'Cmu1  ..  _  cuMgR-t-  M 3i 

da  (i/2)pu2(2b)2 


r 


In  the  ultimate  formulation  to  be  obtained  herein,  it  will  be  useful  to 
express  the  aerodynamic  damping  parameter  in  terms  of  the  work  per  cycle 
of  motion,  since  the  latter  is  a  quantity  which  is  most  easily  measured 
from  the  available  test  data.  To  accomplish  this,  a  few  manipulations  are 
necessary.  First,  Eq.  (7b)  of  Appendix  I  will  be  rewritten  in  the  form 


CMut  -  ^  1  ^ 


(3b) 


where 


/iT(k,Q)  =  [moj  -  (Lax  +  +Q)  +  Lh,  (y  +  of] 


(37) 


is  a  function  of  |<  and  a  only.  Now,  substitute  Eq.  (36)  into 
Eq.  (18)  to  obtain 

Cw  =  TT  flyd2 

and  next  substitute  Eq.  (36)  into  the  derivative  of  Eq.  (35)  with  the 
result 

£7  -  *Cmu1 

~q2  '  da 


=  -M 


(36) 


(39) 


Finally,  after  solving  Eq.  (38)  for  fj.y  ,  substitution  into  Eq.  (39) 
yields  the  useful  formula 


Cw 


*->a. 


va 


n  2 


(kO) 


Ultimately,  these  two-dimensional  data  must  be  used  in  the  stability  eval¬ 
uation  of  the  three-dimensional  rotor  system;  hence,  they  must  be  converted 
to  three-dimensional  aerodynamic  damping  form.  This  has  been  done  in 
Appendix  II . 


=  L  =a2  v*  ( i?)  fu(^)  dr] 


(hi) 


where  I't7?)  is  the  spanwise  velocity  ratio,  f a is  the  spanwise 
mode  shape  distribution,  and  t)  is  the  dimensionless  spanwise  variable. 
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SPECIALIZATION  TO  25%  CHORD  PIVOT.  AXIS 


r .  - .°' ^lisr 

Toil  pivoted  about  its  2^  a'°^  *n*  near  the  25?  chord 

will  also  have  its  effective  pivot  axis  at  m  very  be  useful, 

(assuming  a  conventionally  dcbigue  sect^ons  for  the  special  case 

then,  to  convert  the  results  .  *  +  - a.,  nation  of  the  need  for  these 

of  25?  chord  pivot  axis  location,  m  anticipation  oi  the  neea 

specialised  formulas  in  subsequent  sections. 
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Appendix^I  Zr^it  t<*  «■..  MIN*  P«>'  ^  «•  — «•* 
ficient  amplitude  becomes 

-  .  rtf.  Un  a  l»*2) 

CM(jI--2  MaiQ 


and  substitution  into  Eq.  (18)  yields 


TT2 K2_  u_,a2 


the  work  coefficient.  Similarly,  Eq.  Uo)  becomes 


for  the  two-dimensional  damping  coefficient.  In  the  case  ot^  an  P 

IbU.  potential  fl»  past  the  airfoil ,  tne  i^inar^part^f  S,.  W)  »U 
be  substituted  from  Appendix  I  lor  Mai  ^ 

become 


-  — jr-  7rk  a 


Cw  =  -  —  TT^k  a?- 


Ea.  (46)  shows  that  l'or  these  ideal  conditions,  the  work  coefficient 
will  always  he  negative  (i.e.,  the  system  will  always  he  stable),  and  is 
directly  dependent  on  both  k  and  a  .  Similarly,  Eq.  (47)  shows  that 
the  aerodynamic  damping  is  always  positive  and  is  independent  of  amplitude. 
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FINAL  DATA  REDUCTION 


In  the  course  of  tlie  original  test  program,  approx imately  550  test  points 
were  taken;  for  each  of  these  a  Cm- a  loop  was  generated.  Figures  9,  10, 
and  11  contain  fourteen  separate  examples  of  these  Iccps.  In  principle, 

Sq.  (14)  was  used  to  evaluate  the  work  coefficient  for  each  of  these  loops; 
however,  a  planimeter  was  used  to  obtain  the  actual  area.  Each  loop  was 
planimeterea  at  least  three  times,  and  an  average  of  these  three  measure¬ 
ments  was  used  to  represent  the  area  of  a  given  loop.  (The  maximum  error 
incurred  by  this  averaging  procedure  is  estimated  to  be  less  than  3/5). 
Appropriate  conversion  factors  were  used  to  transform  the  raw  area  data 
into  work  coefficient  form,  and  Eq .  (40)  was  then  used  to  convert  these 
work  coefficients  into  two-dimensional  damping  parameter  form. 

The  results  of  these  calculations  are  presented  in  Tables  I  through  III. 

The  gaps  in  these  tables  denote  values  of  frequency  and/or  mean  incidence 
angle  at  which  either  no  data  were  taken  or  no  data  were  reduced.  For 
example,  only  a  limited  amount  of  very  low  frequency  data  was  taken,  and 
these  points  were  restricted  to  mean  incidence  angles  of  aM  =  12  degrees 
and  greater.  Also,  as  Mach  number  was  increased,  some  of  the  higher  fre¬ 
quency  and  higher  incidence  angle  runs  were  omitted.  Occasionally,  a 
double  entry  was  made  in  the  tables  for  a  particular  combination  of  f  and 
dM  .  This  occurred  whenever  a  Cm- a  loop  failed  to  close  on  itself. 

For  example,  in  Table  1,  for  5=6  degrees  at  the  point  aM  =  3  degrees 
and  f  =4  cps  (  k  =  0.1125)  there  are  two  values.  The  loop  for  this 
case  does  not  return  along  the  same  path  after  one  cycle  of  motion.  In 
this  case,  as  in  other  similar  cases,  an  estimate  was  made  of  the  two  most 
likely  loops,  the  areas  of  both  were  taken,  and  both  numbers  were  entered 
in  the  appropriate  table. 

A  representative  group  of  these  reduced  damping  parameter  values  is  shown 
in  Figures  12  and  13.  These  results  were  plotted  directly  from  the  second 
subtable  of  Table  1.  This  pair  of  figures  shows  the  variation  of  Haz 
with  reduced  frequency  k  for  the  entire  range  of  mean  incidence  angles. 
Included  in  each  plot  is  the  predicted  theoretical  variation  of  Ba2  with 
k  ,  taken  from  Eq.  (47)  for  25%  chord  pivot.  It  is  seen  from  Figure  12 
that  the  experimental  values  are  in  good  agreement  with  theox-y  for  Qm  =  0, 
3,  6,  and  9  degrees.  However,  the  data  show  considerable  departure  from 
the  theory  for  incidence  angles  of  0  m  =  12  degrees  and  above;  in  fact, 
there  is  a  range  of  incidence  angles,  12  <  qm  <  21  degrees,  in  which 

the  aerodynamic  damping  is  negative  (indicating  an  instability)  over  some 
portion  of  the  k  range.  Finally,  at  extremely  large  values  r'f  incidence 
angle,  Qm  >  24  degrees,  the  values  of  "  a2  become  positive  once  again 
and  seem  to  approach  the  theoretical  prediction  at  the  highest  incidence 
angles  shown  here. 


There  are  innumerable  ways  of  presenting  the  entire  mass  of  data  contained 
in  Tables  I  through  III.  One  obvious  method  is  to  go  through  each  sub¬ 
table,  plotting  the  results  in  an  analogous  fashion  to  the  results  shown 
in  Figuies  12  and  13,  with  each  set  representing  a  specific  combination  of 
Mach  number  and  torsional  amplitude.  This  procedure  is  unwieldy  in  the 
ultimate  application  of  the  data,  which  is  to  determine  the  stability 
characteristics  of  helicopter  rotor  blades.  Therefore,  ar-  effort  has  been 
made  to  reduce  the  number  of  independent  parameters  by  combining  these 
data  sets  in  various  ways. 

The  method  to  reduce  the  number  of  independent  parameters  which  was  chosen 
for  this  program  utilized  the  grouping  of  data  as  presented  in  Tables  I 
through  III.  Specifically,  all  of  the  data  for  M  =  0.2  were  combined,  and 
an  effort  was  made  to  eliminate  q  as  a  parameter.  Similarly,  the  data 
for  M  =  0.3  and  for  M  =  O.h  were  also  treated  separately.  The  elimination 
of  5  as  a  parameter  was  effected  by  plotting  the  variation  of  H02 
with  h  for  each  mean  incidence  angle  and  by  including  data  for  all  three 
values  of  the  amplitude.  This  decision  to  eliminate  a  base  in  part 

on  the  fact  that  the  theoretical  result  presented  in  Eq.  ( U7 )  1  j2 

is  independent  of  a 

After  all  of  the  data  had  been  plotted  in  this  fashion  (  ^a2  versus 

k  for  each  aM  ,  for  each  M  ,  with  the  data  for  ail  values  of  a  on 

each  plot),  a  set  of  curves  was  faired  through  the  data  for  each  set  of 

parameter  values.  A  tabulation  of  the  coordinates  of  the  faired  curves 

was  then  made  for  a  large  n'linber  of  conveniently  spaced  values  of  K  at 

the  given  value3  of  aM  ,  ana  a  set  of  cross  plots  of  2  versus 

was  made  using  these  tabulated  values.  Once  again,  a  set  of  curves 
was  faired  through  the  data,  and  the  coordinates  were  taoulated.  this  time 
for  a  closer  spacing  of  incidence  angles  than  were  previously  available 
(every  degree  instead  of  every  3  degrees).  This  process  was  repeated 
sufficiently  often  to  reduce  the  scatter  of  the  coordinates  of  the  faired 
curves  to  a  negligible  amount;  in  effect,  a  graphical  iteration  process 
was  employed  to  smooth  the  data.  (The  reason  for  adopting  this  procedure 
was  to  provide  a  smoothed  set  of  unsteady  aerodynamic  data  suitable  for 
bivariant  interpolation  on  a  high-speed  digital  computer.}  The  results  of 
this  pi'oeess  are  presented  in  Tables  IV ,  V ,  and  VI  for  Mach  numbers  M  - 
0.2,  0.3.  and  O.h,  respectively. 

A  comparison  of  the  original  data  and  the  final  smoothed  data  is  shown  in 
Figures  lU,  15,  and  16.  These  figures  contain  almost  all  of  the  aerody¬ 
namic  damping  data  for  M  ~  0.2  which  were  originally  obtained  from  the 
Cm  ~  a  loops  and  were  presented  in  Table  1.  The  faired  curves  were  taken 
from  Table  IV.  Similar  plots  could  be  constructed  for  M  -  0.3  and  O.h 
from  the  original  data  in  Tables  II  arid  111  and  the  smoothed  data  in  Tables 
V  and  VI.  However,  in  the  interest,  of  brevity,  they  are  not  included  here. 

A  number  of  conclusions  may  be  reached  by  studying  the  results  presented 
In  these  figures,  as  enumerated  below.  It  is  seen  from  Figure  1^  that  the 
experimental  results  are  generally  in  good  agreement  with  the  theory  at  low 
mean  incidence  angle  (  S  9  deg)  for  all  amplitudes.  As  the  mean  inci¬ 
dence  angle  is  increased  to  values  greater  than  aM  =  12  degrees,  the 


aerodynamic  damping  is  seen  to  depart  from  the  theory  and,  in  some  cases, 
t.n  become  negative  over  some  range  of  k  .  A  study  of  Figures  lk  and  .15 
i  s quite  instructive  in  explaining  the  apparent  contradiction  noted  earlier 
in  the  behavior  of  with  increasing  k  .  In  these  figures,  it  is 

seen  that  for  sufficiently  large  k  ,  an  initially  unstable  or  marginally 
stable  condition  will  become  more  stable  as  k  is  increased,  a.s  stated 
earlier  in  this  reper..  However,  at  very  low  values  of  k  ,  the  aerody¬ 
namic  damping  initially  increases  from  a  zero  value  at  k  =  C  and  reaches 
a  small  positive  value  before  becoming  negative  at  a  higher  value  of  K 
This  vac  pointed  out  earlier  in  connection  with  the  behavior  of  the  loops 
in  Figure  11.  Finally,  as  shown  in  Figure  15  for  sufficiently  large  values 
of  Cl m  ,  the  system  once  again  exhibits  a  stable  behavior. 

A  similar  behavior  in  the  variation  of  « with  Qy  is  seen  to  exist 
for  intermediate  values  of  k  (Figure  16).  This  figure  shows  the  aerody¬ 
namic  damping  to  be  initially  stable  at  low  incidence  angle,  unstable  in 
the  approximate  range  cf  12  degrees  <  <2m  <  22  degrees,  and  then  stable 

once  again.  Thus,  then  appears  to  be  a  "pocket"  of  instability  surrounded 
by  stable  regions  at  botn  higher  and  lower  incidence  angles,  and  at  lower 
values  uf  reduced  frequency. 

This  isolated  region  of  instability  is  graphically  illustrated  in  Figures 
17,  lti,  and  19.  In  these  figures,  the  aerodynamic  damping  is  represented 
as  a  surface  relative  to  the  k  ,  -plane  for  Mach  numbers  of  ii  = 

0.2,  0.3,  and  0.1*.  These  aerodynamic  damping  surfaces  are  presented  in 
both  an  isometric  view  and  a  near  platform  view  in  each  figure.  The  super¬ 
imposed  rent angular'  grid  lines  on  each  surface  represent  constant  values 
of  either  k  or  aM  ,  and  the  dark-colored  region  lying  below  the  k  , 
-plane  is  a  region  of  negative  aerodynamic  damping  and,  hence,  repre¬ 
sents  11  potentially  unstable  aerodynamic  condition. 

CONVERSION  OF  REFERENCE  1  DATA 


A  number  of  years  ago,  an  extensive  experimental  program  was  carried  out 
as  reported  in  Reference  3  to  determine  the  Variations  in  lift  and  moment 
coefficient  on  a  12%  thick  airfoil  oscillating  in  both  pitch  and  transla¬ 
tion  over  a  wide  range  of  incidence  angles.  This  work  was  limited  in  some 
respects  by  the  restriction  to  low  free-stream  velocities  and  by  the  fact 
that  the  pivot  axis  for  pitching  motions  was  located  at  the  375!  chord  sta¬ 
tion  (  a  -  -  0.26).  'However,  the  type  of  airfoil  employed  in  the  tests 
and  the  range  of  parameters  over  which  experimental  data  were  obtained  make 
the  Reference  1  results  potentially  quite  useful  in  the  present  analysis. 
The  remainder  of  this  section  is  devoted  tc  the  transformation  necessary 
to  convert  these  coefficients  to  «  form  which  can  be  compared  directly  with 
the  results  of  the  present  analysis. 

The  experimental  date,  tabulated  i r >  Reference  I  are  presented  in  the  form  of 
an  ampj i nude  and  prase  angle  for  each  component .  Specifically,  the  quan¬ 
tities  involved  are  CLT  ,  <£>LT  ,  cUP  ,  4>  Lp  »  CMt  ,  <PMJ  >  CMP  ,  4>MP  , 
where  C  represents  the  aropli tv.de  of  the  coefficient;  <£  ,  the  phase 
angle  by  which  the  force  or  moment,  leads  the  motion;  and  the  subscripts 
LT  ,  lp  ,  n/,i  ,  MP  ,  the  lift  due  to  translation,  lift  due  to  pitch, 


moment  due  to  translation ,  and  moment  due  to  pitch,  respectively.  It  will 
be  shown  below  that  all  four  coefficients  and  their  phase  angles  are  needed 
to  convert  the  Reference  1  data  to  a  form  comparable  with  the  data  pre¬ 
sented  in  this  volume.  This  is  because  the  results  given  in  Reference  1 
were  obtained  for  a  37 ^  chord  pivot  axis,  whereas  the  present  repor"  is 
concerned  with  a  pivot  axis  at  the  25?  chord  station.  Other  major  dif¬ 
ferences  between  the  Reference  1  and  the  present  results  involve  the  def¬ 
initions  of  the  coefficients  and  the  fact  tnp-t  the  measured  unsteady  coef¬ 
ficients  of  Reference  1  were  not  independent  of  t ne  amplitudes  of  oscilla¬ 
tion.  These  differences  will  be  resolved  in  the  course  of  the  analysis. 

In  Reference  1  the  lift  and  moment  coefficients  are  defined  by  the  equa¬ 
tions 


— V- 

2?uzb 


(45; 


where,  for  example, 

clt"  r  ClT(cos0LT  t  i  sm<*>LT) 


(50) 


and  where  the  superscript  H  in  Eqs.  (48)  and  (49)  is  used  to  denote  the 
coefficient  form  employed  in  Reference  1.  Eq.  (48)  and  (49)  may  be  solved 
for  the  unsteady  lift  and  moment  as 

ZJ  --  2/>U2b  (CLT*+  cLP*)  (51) 


Mu  '  2.P  U2b2  ( CMT  +CMP  )  (  52) 


Note  that  neither  the  dimensionless  bending  deflation  amplitude  Fi  -  h/b  nor 
the  twist  deflection  a  appears  on  the  right  side  of  these  equations. 

As  stated  earlier,  this  is  because  the  Reference  1  coefficients  are  not  in¬ 
dependent  of  the  amplitudes  of  motion.  The  actual  values  of  the  parameters 
used  in  Reference  1  will  be  substituted  into  these  equations  later  in  the 
development . 


0 


The  unsteady  lift  and  moment  ampj ituo.es  associated  with  both  translation 
and  pitch  about  a  pivot  axis  at  q  may  be  taken  from  References  8  of  13 . 

Lu*  =  ^bV{Lhh  +[l0 -(f -0)Lh]d}  (53) 

fj[*  =  TTpb^2  {[Mh  -  (^' +  o)Lh]h 

+  [m0-(  ^-  +  Q){l_0  +  Mh)  +  (7+a)Zl.fijo} 

(5>M 

A  combination  of  £qe .  ( 51)  and  (53)  and  bqs .  (52)  and  ( )  yields  the 
relationships 

Lh  h  +  I  LU  -  (  j  +  o)  Lh]  5  :  f  Clt  +  C lp  )  (55) 

and  ■'  '  1 

[Mh_(^  +  a)>-h]b  +  [^q  ~{  2  +  o)(La  +  Mh)  +  (sy  +  a)  i_h  ja 

=  (cmt  +  Cm*  )  (56) 

After  separating  the  motions  into  single-degree-of- freedom  oscillations  in 
bending  and  torsion,  Eqs.  (55)  and  (5 6)  may  be  recast  in  an  equivalent  form 
as  four  equations,  as  follows: 
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[mq-(1  +  q) 

(lq  +  Mh)  +  ( 

+  a)  Lh]  a  =  cmp 

(CO) 
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Tissue,  in  turn,  may  be  solved  for  the  standard  coefficients  ,  Lq 

M  h  »  M  (j  : 


The  quantities  on  the  right  side  of  Eqs .  (61)  through  (64)  represent  the 
experimental  values  obtained  in  Reference  i  .  Numerical  values  must  be  in¬ 
serted  for  the  pivot  axis  location,  a  =  -  0.26,  and  for  the  dimensionless 
displacements,  b  =  h/b  “  0-9/5-808  =  0.155,  and  a  =  6.08  7r/l80  =  0.106. 
The  quantities  found  on  the  left  sides  of  these  equations  will  also  repre¬ 
sent  the  experimental  values;  however,  these  quantities  will  be  independent 
of  pivot  axis  location,  at  least  according  to  the  linear  theory  employed 
here. 

To  use  the  Reference  1  data  in  the  present  analysis,  the  imaginary  part  of 
Eq.  (64)  must  be  substituted  into  Eq.  (44).  After  inserting  the  value 
0  =  -  0.26,  the  result  for  the  two-dimensional  aerodynamic  damping  is 

Sa  =  ~  Tt[cmp  sin*MP  +  °.24CLPsin<£  J 

“  ■2}pL[CMT  sin^MT  +  0  24 CLT  sin<pLTj  (65) 


This  equation  was  used  to  evaluate  _  ~a2  from  the  Reference  1  coef¬ 
ficients  after  inserting  the  values  h  =  0.155  and  a  =  0.106. 


COMPARISON  OF  PRESENT  RESULTS  WITH  REFERENCE  1  DATA 

Before  these  results  could  be  applied  in  the  present  stability  analysis, 
some  basis  of  comparison  had  to  be  established  for  the  present  data  and  the 
Reference  1  data,  as  described  bflow.  Three  configurations  vere  tested,  in 
Reference  1  — the  blunt  wing,  the  intermediate  wing,  and  the  sharp  wing — 
in  which  the  profile:'  differed  only  in  the  vicinity  of  the  leading  edge. 
Since  there  is  a  slight  difference  in  both  the  steady-state  (i.e.,  stalling 
angle)  and  nonsteadv  behavior  of  these  configurations,  the  data  for  each 
are  presented  separately  in  Reference  1  as  functions  of  mean  incidence 
angle,  ,  and  reduced  frequency,  k  .  It  was  decided  to  correlate 

these  results  with  one  another  and  with  the  present  data  by  means  of  a 
stall  angle  parameter, 

°Y  =  aM/QS  (66> 


The  angle  as  is  the  stalling  angle  for  the  particular  configuration  under 
consideration.  This  varied  from  configuration  to  cojif 5 guration  in  Reference 
1,  in  the  present  case  it  was  taken  to  be  Qj  =  13  degrees,  which  is  the 
incidence  angle  at  which  the  slope  of  the  steady-state  normal  force  curve 
becomes  horizontal  in  Figure  8.  This  choice  of  a  correlation  parameter 
appeared  to  be  a  logical  one  in  view  of  the  strong  effect  the  stall  point 
has  on  unsteady  response.  It  can  be  seen  from  Figure  16,  for  example,  that 
the  atrodynnc.ic  dumping  is  nearly  constant,  below  the  stalling  angle  and 
begins  to  vary  only  when  the  combination  of  mean  incidence  angle  and  ampli¬ 
tude  of  motion  cause  the  airfoil  to  be  affected  by  stalling  effects. 

An  in  the  case  of  the  present  results,  the  Reference  1  data  were  bivariantly 
dependent  on  both  Q  and  k  ;  a  graphical  iteration  scheme,  similar  to  the 
method  used  on  the  present  data,  was  employed  to  smooth  these  data.  Equa¬ 
tion  (65)  was  then  used  to  calculate  the  two-dimensional  aerodynamic  damping 
as  a  function,  this  time,  of  stall  angle  parameter,  crT  t  and  reduced  fre¬ 
quency,  k  .  These  results  are  tabulated  in  Table  VII. 

A  comparison  of  the  two-dimensional  aerodynamic  damping  results  from  the 
present  data  and  from  the  Reference  1  data  is  presented  in  Figures  20  and 
21.  It  is  seen  from  these  figures  that  the  agreement  between  the  two  sets 
of  data  is  quite  good  at  row  values  of  k  ,  is  excellent  at  intermediate 
values  of  k  ,  and  is  qualitatively  good  at  the  highest  value  of  k  .  Of 
particular  significance  is  the  fact  that  both  sets  of  data  predict  essen¬ 
tially  the  same  two-dimensional  stability  boundary  at  moderate  values  of 
k  in  the  range  0,2  £  k  <  0-3  and,  in  addition,  predict  the  same  damping 
level  as  a  function  of  .  Thus,  it  is  seen  that  the  pivot  axis  trans¬ 

formation  procedure  is  valid,  and  the  use  of  the  stall  angle  parameter  as 
the  basis  for  correlation  is  also  valid,  at  least  within  the  variable  ranges 
considered  herein. 


One  further  comparison  of  the  two-dimensional  aerodynamic  damping  parameter 


ccent  stuny  and  those  of  Reference  it. 

The  aerodynamic  damping  curve  for  25?  chord  pivot  and  U/ buj =  3.0  (i.e., 
for  k  =  0.333)  was  taken  from  Figure  85  of  Reference  lU  (which,  in  turn, 
was  obtained  from  a  cross  plot  of  previous  results  presented  in  Reference 
12)  and  was  compared  with  the  curve  for  M  =  0.2,  k  -  0-3375  from  Figure 
16  of  this  volume.  A  direct  comparison  was  impossible  to  make  since  the 
airfoil  section  of  Reference  12  was  a  thin  (less  than  O.Oii  thickness  ratio) 
MACA  65-series  profile  which  undoubtedly  did  not  have  the  same  stalling 
characteristics  as  the  NACA  0012  airfoil  used  herein.  Unfortunately,  the 
static  stall  angle  of  this  65-series  airfoil  was  not  provided  in  Reference 
12.  Hence,  to  provide  some  sort  of  comparison  of  the  Reference  lU  results 
and  the  present  results,  it  was  necessary  to  improvise  a  basis  of  correla¬ 
tion.  To  this  end,  it  was  decided  to  scale  the  abscissa  in  Figure  lH  of 
Reference  lU  to  force  a  coincidence  of  the  first  zero  crossing  of  Sa2 
for  the  25?  chord  pivot  curve  with  the  first  zero  crossing  of  the  curve  in 
Figure  16.  This  comparison  is  presented  in  Figure  22,  in  which  the  abscissa 
has  been  converted  to  stall  angle  parameter  form. 


It  is  obvious  that  this  comparison  can  have  only  qualitative  value,  in 
view  of  the  liberties  taken  with  the  scale  of  the  abscissa  values  from 
Reference  lU,  nevertheless,  the  comparison  shows  a  general  agreement 
between  the  two  results ,  including  the  extent  of  the  negative  damping 
region  and  the  recovery  at  higher  values  of  oy  ,  despite  the  large  dif¬ 
ferences  in  the  airfoil  configurations  and  the  test  procedures  employed. 


<j 


ROTOR  STABILITY  CALCULATIONS 
DESCRIPTION  OF  TYPICAL  ROTOR  CONFIGURATION 

The  Sikorsky  S-61F  rotor  was  chosen  for  the  application  of  the  stability 
analysis  described  in  previous  sections  of  this  report.  The  rotor  blade 
had  a  total  chord  of  1.921  feet  and  extended  from  a  root  radius  of  ro  = 
8,k38  feet  to  a  tip  radius  of  rT  =  31.0  feet.  Performance  data  for  this 
configuration  were  calculated  for  four  flight  conditions:  forward  speed  of 
16.1  knots  and  10,660  pounds  gross  weight,  165  knots  and  12,k60  pounds,  lb9 
knots  and  18,200  pounds,  and  210  knots  and  16,200  pounds.  The  tabulated 
data  provided  for  each  condition  consisted  of  radial  and  azimuthal  varia¬ 
tions  in  incidence  angle  and  Mach  number. 

It  was  assumed  that  the  rotor  blade  was  capable  of  responding  to  an  in 
nitesimal  disturbance  in  its  fundamental  torsional  mode  at  every  azimuth, 
position,  and  the  stability  analysis  was  an  ployed  to  determine  the  damping 
of  this  incipient  vibration.  The  fundamental  torsional  frequency  of  the 
S-61F  rotor  blade  is  f  =  27-3  cps  or  u>  ~  171.^  rad/sec;  this  quantity, 
together  with  the  semichord  dimension  b  =  0.76  feet,  was  used  in  Eq.  (7) 
to  calculate  the  numerator  of  the  reduced  frequency  parameter ,  k  .At 
each  azimuthal  station  selected  for  the  analysis,  the  radial  variation  in 
Mach  number  was  converted  to  a  comparable  radial  variation  in  velocity  by 
multiplying  each  value  by  the  assumed  local  speed  of  sound,  1100  ft/sec. 

This  was  used  to  calculate  the  radial  variation  in  k  which,  together  with 
the  given  radial  variation  in  a  ,  was  used  to  interpolate  the  appropriate 
aerodynamic  damping  table  for  the  required  value  of  Sq2  a't  each  radial 
station.  (Polar  plots  of  typical  k  -  and  a  -  distributions  over  the 
rotor  disc  are  shown  in  Figure  23.)  The  reference  velocity,  Ut  >  was 
chosen  to  be  the  tip  value  of  the  velocity  at  the  zero  azimuth  station, 

4/  =■•■  0°,  and  it  was  used  to  compute  the  velocity  ratio,  2/(p)  ,  defined  by 
Eq.  (98)  in  Appendix  II.  The  computed  torsional  mode  shape  function,  , 

was  calculated  for  the  S-61F  helfcopter  rotor  blade  and  is  presented  In  Table 
VIII.  With  the  introduction  of  these  quantities  into  Eq.  (Ll),  the  span- 
wise  integral  was  calculated  and  the  three-dimensional  damping  parameter, 

Hqj  (vM  ,  was  obtained  for  each  value  of  azimuth  angle. 

STABILITY  ANALYSIS 

Three  different  procedures  were  employed  in  extracting  information  from  the 
tables.  In  the  first  case,  only  the  M  -  0.2  table  was  used  and  the  actual 
Mach  number  was  used  only  to  calculate  U  .  If  the  value  of  either  k 
or  a  exceeded  the  tabulated  value  in  Table  IV  in  a  corner  value  of  the 
damping  parameter  was  used;  i.c.,  the  oversize  parameter  was  artifically 
returned  to  the  maximum  tabulated  value,  and  the  interpolation  with  respect 
to  the  other  parameter  was  performed. 

A  second  case  involved  data  for  all  three  Mach  numbers  from  Tables  IV 
through  VI.  Here  the  actual  Mach  number  was  used  to  interpolate  between 
tables  as  well  as  to  provide  values  of  th  ’  velocity,  U  .  Corner  values 
were  also  taken  here;  however,  because  of  the  truncation  of  the  tables  for 
the  higher  Mach  numbers,  this  procedure  led  to  somewhat,  unsatisfactory 
compromises . 


The  third  case  utilized  the  Reference  1  data  presented  in  Table  VII.  Once 
again,  the  Mach  number  was  used  to  calculate  U  only,  as  in  the  first 
case;  and,  once  again,  corner  values  were  employed. 

Stability  analyses  for  the  four  flight  conditions  mentioned  earlier  were 
carried  out  by  using  each  of  these  three  options.  The  results  are  pre¬ 
sented  in  Figures  24  and  25,  in  which  the  variation  of  the  three-dimensional 
aerodynamic  damping  parameter,  a i  ,  with  azimuth  angle,  ^  ,  is  pre¬ 

sented.  In  each  figure  the  abscissa  has  been  extended  by  It 0  degrees  beyond 
a  full  cycle  to  show  more  clearly  the  behavior  of  the  damping  in  the  region 
surrounding  =0  degrees. 

In  Figure  24,  for  a  flight  speed  of  165  knots  and  moderate  gross  weights, 
all  three  sets  of  data  are  in  good  agreement  and  predict  a  stable  operation 
over  all  values  of  ^  .  The  results  from  the  present  data  are  somewhat 

more  conservative  than  the  Reference  1  results  for  the  advancing  blade 
region  (  0°  -  ^  -  180°) ,  and  all  three  inputs  are  in  very  close  agreement 
over  the  retreating  blade  region  (l8o°  $  'I'  £  360°).  As  expected,  the 
damping  is  greater  in  magnitude  for  0°  -  ^  i  l8o°  than  it  is  for 
l80°  <<p  -  360°.  This  is  because  the  dynamic  pressure  attains  a  maximum 

value  in  this  region;  hence,  the  square  of  the  velocity  ratio,  which  weights 
the  integrand  in  Eq.  (h±)  will  also  attain  a  maximum  value. 

In  the  upper  portion  of  Figure  25  the  flight  speed  has  been  reduced  slight¬ 
ly,  but  the  gross  weight  has  experienced  a  significant  increase,  relative 
to  the  previous  two  figures.  Once  again,  the  results  of  the  three  damping- 
table  options  are  in  generally  good  agreement;  and,  once  again,  the  Ref¬ 
erence  1  data  yield  less  conservative  results  than  the  present  data,  in  the 
range  0°  i  t  -  180°.  It  is  seen  that  in  the  region  of  the  retreating 
blade,  the  system  damping  decreases  to  very  small  values,  which  implies 
reduced  flutter  margin  for  the  retreating  blade. 

Finally,  in  the  lower  portion  of  Figure  25,  both  flight  speed  and  gross 
weight  are  considerably  increased  relative  to  the  previous  three  cases.  As 
shown  in  Figure  23,  an  extensive  region  of  the  rotor  disc  is  simultaneously 
operating  at  large  values  of  k  and  at  incidence  angles  exceeding  the 
steady-state  stall  angle  (  Qs  =13  degrees).  As  a  result,  both  the  present 
data  for  M  =  0.2  only  and  the  Reference  1  data  predict  a  region  of  large 
negative  damping  over  the  approximate  range  340°  S  \j/  <  366°.  This  be¬ 
havior  will  be  considered  in  detail  in  a  later  section.  (Note  that  the 
present  data  for  M  =  0.2  only  also  predict  a  shallow  region  of  negative 
damping  over  the  approximate  range  233°  -  -  251°.  This,  too,  will  be 

studied  in  more  detail  later.)  First,  it  should  be  noted  that  although  the 
present  data  taken  from  all  Mach  numbers  do  not  agree  with  the  other  two 
sources  in  predicting  the  instability  for  the  retreating  blade,  they  do  agree 
elsewhere  for  the  advancing  blade.  This  disagreement  is  caused  in  part  by 
the  method  adopted  in  the  computer  program  in  choosing  the  corner  values. 

The  corner  values  used  in  the  case  of  free  interpolation  among  Mach  numbers 
are  determined  by  the  tabular  value  of  M  nearest  to  the  actual  local  value 
of  M.  It  is  seen  from  Figures  18  and  19  that  available  damping  information 
is  truncated  for  both  a  and  k  as  M  is  increased.  Hence,  although  k 
and  a  may  have  large  values  on  the  retreating  blade,  if  the  Mach  number 


is  high  enough,  the  corner  value  imposed  by  the  interpolation  scheme  may 
yic'd  a  value  of  the  damping  which  is  considerably  removed  from  the  actual 
damping  for  the  given  conditions.  It  is  felt  that  this  compromise  limits 
the  usefulness  of  the  free  Mach  number  interpolation  scheme;  and  the  re¬ 
mainder  of  the  resuJ ts  considered  herein  will  be  restricted  to  the  use  of 
the  present  data  for  M  =  0.2  only,  with  occasional  reference  to  the  use  of 
the  Reference  1  data. 

The  actual  effect  of  there  potentially  unstable  regions  on  the  stall  flut¬ 
ter  response  of  the  rotor  is  probably  not  too  serious  for  the  few  cases 
considered  so  far  in  this  report.  This  is  because  1  he  torsional  vibration 
frequency  of  27-3  ops  corresponds  to  an  3-eyclo-per-revolution  torsional 
motion;  hence,  each  complete  cycle  will  extend  over  ^5  degrees  of  azimuth. 
The  largest  of  the  two  unstable  regions  predicted  by  the  use  of  the  present 
data  l'or  M  =  0.2  extends  for  only  25  degrees  of  azimuth;  therefore  it  is 
incapable  of  exciting  more  than  a  hall'  cycle  of  torsional  response.  There¬ 
after,  the  aerodynamic  damping  returns  to  a  large  positive  value,  which  is 
sufficient  to  damp  out  the  motion. 

COMPARISON  WITH  RESULTS  OP  REFERENCE  2 


After  the  work  reported  herein  vi.i  v-1'  under  way,  it  was  found  that  a 
parallel  study  had  been  carried  out.  Reference  'u  which  both  an  oscill¬ 
atory  limit-cycle  experiment  was  conducted  and  a  stability  analysis  based 
on  available  aerodynamic  data  was  carried  out.  The  significant  result  of 
Reference  2  is  the  fact  that  a  region  of  instability  can  exis ;■  for  the 
retreating  blade.  This  is  shown  in  Figure  13  of  Reference  2  and  it  is 
interesting  to  note  that  the  damping  variation  with  p  portrayed  in  this 
figure  first  increases  with  tp  for  the  advancing  blade  and  red' he'  a 
maximum  value  at  approximately  V  =  00  degrees,  after  which  it  creases 
to  zero  at  approximately  y  =  225  degrees ,  which  indicates  an  instability 
over  the  approximate  range  225° <  4  <3C0°.  '  is  bears  a  •  try  strong  re¬ 

semblance  to  the  results  shown  in  the  lower  portion  of  ?  Lgiri  5  of  the 
present  volume,  although  there  are  certain  significant  diffe  jees  between 
the  two  results,  particularly  in  the  region  of  the  retreating  blade. 
Specifically,  the  present  i-tjsul  i  indicates  two  narrow,  unstab.'  s  regions 
separated  by  a  stable  region,  while  the  Reference  2  result  pr<  lets  a 
rather  broad,  continuous  region  of  instability.  (Note  that  the  result 
based  on  Reference  1  data  in  Figure  25  also  agrees  closely  with  the  result 
from  the  present  data;  the  deep  instability  over  the  range  3^0°  <  <p  < 

366°  is  matched  almost  exactly, and  there  is  a  tendency  toward  instability 
in  the  vicinity  of  ip  =  2)i0  degrees.) 

The  general  agreement  between  the  present  results  and  tho.  e  of  Reference  2 
is  very  encouraging  in  view  of  the  f at t  that,  although  these  were  parallel 
studies,  they  were  carried  out  independently  of  one  another  and  utilized 
distinctly  different  sets  of  input  data,  Leith  i\.,r  i  otor  performance  charac- 
teristics  and  for  unsteady  aerodynamic  damping  data. 

The  differences  in  input  data  lead  directly  to  the  differences  in  the  sta¬ 
bility  results  on  the  retreating  blase.  First,  the  damping  data  used  in 
Reference  2  were  synthesized  iron,  earlier  results  in  Reference  li ,  12,  and 


15-  Although  these  data  were  corrected  J'or  static  stalling  angles  anti 
rotational  axis  location,  the  basic  profile  shapes  were  undoubtedly  d'f- 
ferent  l’rom  the  MCA  0012  profile  used  herein  and  could  conceivably  lead  to 
differences  in  tlie  coefficient  values.  In  particular,  the  first  two-dimen¬ 
sional  damping  curve  to  become  negative  in  Figure  8  of  Reference  2  is  that 
for  o—8  degrees;  whereas  in  the  present  study,  the  first  to  become  neg¬ 
ative  is  that  for  o  =13  degrees,  as  seen  in  either  Figure  lb  or  in  Table 
IV.  These  differences  in  damping  coefficient  behavior  are  certainly  re¬ 
flected  in  the  three-dimensional  damping  response  differences  noted  above. 

A  second  important  difference  between  the  two  results  lies  in  the  rotor 
loading  characteristics  considered  in  each  case.  In  Reference  2,  the  max¬ 
imum  blade  incidence  angle  is  o  =  16  degrees  on  the  retreating  blade, 
which  yields  almost  the  maximum  possible  value  of  negative  aerodynamic 
damping.  In  contrast  to  this,  the  maximum  rotor  incidence  angle  for  the 
present  case  is  a  =  bO  degrees  on  the  retreating  blade  (see  Figure  23) 
which  is  considerably  beyond  the  point  of  maximum  negative  damping  and,  in 
fact,  is  also  considerably  beyond  the  point  where  the  two-dimensional  damp¬ 
ing  becomes  positive  or.ee  again.  Thus,  the  double  instability  predicted  in 
Figure  25  is  explained  quite  simply  by  a  careful  study  of  Figure  23.  In 
particular.  Figure  ?3  indicates  that  an  initial  region  of  instability  may 
'  >  ;ountered  in  the  third  quadrant  as  the  incidence  angle  increases 
thruuj.t  the  negative  damping  regime,  followed  by  a  stable  region  as  the 
positive  damping  regime  is  penetrated  at  high  incidence  angles.  A  second 
region  of  instability  may  be  encountered  near  the  interface  of  the  fourth 
and  first  quadrants,  as  the  incidence  angle  decreases  once  agiin  through 
the  negative  damping  regime. 

STABILITY  IMPLICATIONS  OF  LOADING  VARIATIONS 

The  observations  of  botli  the  similarities  and  differences  between  the  pres¬ 
ent  work  and  the  results  of  Reference  2  led  rather  naturally  to  the  question 
of  rotor  loading  variations  and  to  what  their  effect  on  rotor  stability 
might  be.  t ;  eifioally,  it  was  decided  to  reduce  the  effective  loading  of 
the  maximum  .uat.  condition  of  the  S-tilF  to  determine  whether  or  not  the 
system  would  experience  the  same  extended  region  of  instability  as  that 
described  in  Reference  2.  Rather  than  employing  an  exact  rotor  performance 
calculation  for  each  new  loading  condition  desired,  it  was  decided  instead 
to  multiply  the  entire  incidence  angle  distribution  over  the  rotor  by  a 
constant  factor.  In  other  words,  a  factor  of  0.5  might  be  applied  to  the 
incidence  angle  distribution  over  the  rotor  disc,  whereupon  all  incidence 
angles  would  be  halved  before  the  table  search  for  damping  values  would  be 
employed.  No  changes  were  made  in  the  k-distribution .  Although  this  pro¬ 
cedure  neglected  the  effect  of  loading  changes  on  rotor  inflow,  it  was  felt 
that  any  error  incurrre'.  would  be  of  second  order  in  relation  to  the  ob¬ 
served  changes  in  stability  characteristics. 

The  values  of  the  incidence  angle  multiplier  chosen  for  thin  study  were 
0.9,  0.8,  0.7,  0.667,  0.5  and  0.375.  The  results  obtained  with  obese  mul¬ 
tipliers  are  presented  in  Figure  26.  It  is  clear,  from  the  bottom  portion 
of  Figure  26  that  the  incidence  angle  reduction  produced  by  the  0.8  and  0.9 
multipliers  was  insufficient  to  remove  the  stable  region  lying  between  she 


unstable  extremes,  although  the  two  unstable  regions  tended  to  move  toward 
one  another.  Furthermore,  the  azimuthal  extent  of  each  unstable  region  in¬ 
creased  with  decreasing  multiplier  level. 

In  the  center  portion,  the  0.667  and  0-7  multipliers  both  yielded  contin¬ 
uously  unstable  regions.  In  both  of  these  cases,  it  appears  that  the  mul- 
tiplier  caused  a  sufficient  reduction  in  Q  to  prevent  the  incidence  angle 
from  exceeding  the  upper  limit  of  the  unstable  region.  Although  the  ini¬ 
tial  and  final  azimuth  angles  for  instability  are  somewhat  less  than  they 
were  for  the  original  data,  the  total  extent  of  the  unstable  region  is  con¬ 
siderably  greater  for  the  0.7  multiplier  than  for  the  original  data.  Spe¬ 
cifically,  it  encompasses  9!?  degrees  of  azimuth,  for  an  8-cycle-per- 
revolution  torsional  motion,  this  amounts  to  slightly  more  than  two  full 
cycles  of  motion  which  can  be  excited  by  this  unstable  region.  This  be¬ 
havior  is  now  quite  consistent  with  the  result  given  in  Reference  2. 

Finally,  in  the  top  portion  of  Figure  26,  the  multipliers  0.5  and  0.375  are 
small  enough  to  reduce  the  incidence  angle  to  levels  below  the  unstable 
limit.  In  fact,  use  of  the  0.375  multiplier  yields  a  near  sinusoidal  re¬ 
sponse,  characteristic  of  potential  flow  behavior.  No  stall  flutter  would 
be  possible  with  these  multipliers. 

A  further  clarification  of  these  stability  characteristics  is  afforded  by 
Figure  27  in  which  the  k -  a  trajectory  for  one  cycle  of  motion  for  each 
of  three  multiplier  values  has  been  superimposed  on  the  top  view  of  the 
aerodynamic  damping  suri'ace  for  M  =  0.2  only.  This  figure  clearly  shows 
that,  in  the  case  of  the  original  data,  the  trajectory  passes  completely 
through  the  unstable  region  twice  and  emerges  into  the  stable  regions  on 
either  side  of  it  each  time.  This  accounts  for  the  double  region  of  in¬ 
stability  for  this  case.  It  is  also  seen  that  the  trajectory  for  a  mul¬ 
tiplier  of  0.7  remains  within  the  unstable  region  over  a  substantial  por¬ 
tion  of  each  cycle,  which  produces  the  extensive  region  of  instability. 
Finally,  the  trajectory  for  a  multiplier  of  0.5  is  seen  to  miss  the  un¬ 
stable  region  completely. 


RESULTS 


With  regard  to  the  two-dinens i onul  damp! ng  parameter  in  pitch  for  the  NACA 

0012  airfoil, the  following  results  were  obtained: 

1.  The  measured  values  of  the  two-dimensional  ueiodynamic  damping  par¬ 
ameter,  Ha?  ,  are  in  good  agreement  with  potential  flow  theory  for 
low  mean  incidence  angles,  Qy  <  90  .  and  der-'  "t  from  tlie  theory  for 
Qm  ^  12°  • 

2.  At  low  incidence  angles,  ”ct2  increases  linearly  wit))  reduced 
1 requency  parameter  ,  k 

3.  At  high  incidence  angles,  ~az  i- 'tially  increases  with  k  , 
then  decreases  to  a  negative  level  (implying  instability),  and  finally 
increases  again. 

U.  Die  Reference  1  and  the  present  data  are  in  generally  good  agi  ement 
over  a  substantial  range  of  both  k  and  a M  .  Both  sets  of  data 
•predict  essentially  the  same  two-dimensional  stability  boundary  at 
moderate  values  of  k  .  Othe)  results  for  much  different  profile 
shapes  are  in  good  qualitative  agreement  with  the  present  data. 

With  regard  to  the  three-diironr.l onal  damning  in  pitch  and  the  attendant 

stability  of  an  S-61F  helicopter  rotor  blade,  the  following  results  were 

obtained: 

5-  At.  low  blade  loading  and  flight  sjoed,  the  rotor  blade  is  stable  over 
the  entire  disc.  The  damping  in  pitch  is  greater  for  the  advancing 
blade  than  for  tne  retreating  blade. 

6.  At  large  blade  loading  and  flight  speed,  the*  rotor  blade  is  still 
stable  over  the  advancing  portion  but  tends  toward  instability  over 
the  retreating  portion  of  the  rotor  disc, 

7-  In  the  case  of  the  largest  loading,  and  flight  speed  considered  thus 
far,  the  predicted  negative  damping  region,  for  the  retreating  blade 
is  confined  to  two  small  pockets,  each  of  which  extends  over  an  in¬ 
sufficient  azimuth  range  to  ex’ite  more  than  one-half  cycle  of  tor¬ 
sional  motion. 

8.  Use.-  of  the  Reference  1  data  yields  slightly  larger  blade  dumping  for 
the  advancing  blade  but  j  reducer  gooo  agreement  with  the  present  data 
for  the  retreat 3  ng  blade,  “ven  in  tne  case  of  r  etreating  Linde  insta- 
bili  -y . 

If  the  loading  is  reduced  (without  chungirig  inflow),  a  condition  can 
he  reached  wherein  two  full  cycles  of  torsional  motion  can  be  excited 
by  the  negative  damping  of  the  retreating  blade. 


9. 


The  results  of  the  present  stability  analysis  are  in  good  qualitative 
agreement  with  the  results  of  Reference  2. 


The  available  two-dimensional  data  for  M  =  ( 
in  the  stability  analysis.  The  data  for  M  1 
they  are  limited  to  the  low-velocity  intompi 


1.3  are  inadequate  for  use 
-0.2  are  sufficient,  but 
~essible  flow  regime. 


TABLE  I.  ORIGINAL  TOO-DIMENSIONAL  AERODYNAMIC 
DAMPING  DATA,  H  ,  FOR  M  =  0.2 


f(cps)  = 
k 

"  0.5 

o.oi4o 

1.0 

0.0281 

2.0 

O.O562 

4.0 

0.1125 

8.0 

0.2250 

12.0 

0-3375 

16.0 

0.4500 

fv  O 

UM 

Pitching 

Amplitude, 

6  =  4° 

0 

0.111 

0.222 

0.386 

0.536 

0.111 

0.261 

0.386 

0.536 

6 

0.092 

0.235 

0  •  373 

0.523 

9 

0.092 

0.261 

0.392 

0.516 

12 

-0.020 

-0,052 

-0.020 

0.072 

0.307 

0.451 

0.431 

15 

-0.020 

-0.183 

-0.235 

-0 . 346 

-.026/-. 065 

.268/.046 

18 

0.536 

-0.686 

0.360 

■785/. 569 

21 

2b 

27 

0 . 477 

0.157 

0.732 

1.124/1.431 

30 

0.307 

0.065 

.523/. 503 

1.131/1.046 

32 

-0.118 

-0.190 

.602/. 91 5 

-.301/-. 810 

aw° 

• 

Pitching  Amplitude, 

5=6° 

0 

0.131 

0.259 

0.^85 

0.636 

3 

.177/.. 119 

0.264 

0.447 

O.619 

6 

0.137 

0.279 

0.462 

0.651 

9 

0.160 

0-714 

0.482 

0 

- 

12 

-0.012 

0  017 

0.029 

-0.044 

0.058 

0.392 

0.6o4 

15 

-0.081 

-0.078 

0.099 

0.055 

-0.177 

.038/. 195 

0.351 

18 

-O.lld 

0.110 

0.01.5 

0.099 

-0 . 360 

-0.674 

-0.293 

21 

1 

0.404 

0.459 

-0.346 

-0-572 

-0.468 

2b 

0.029 

0.020 

O.O58 

0.471 

0.911 

0.813 

0.279 

27 

0.049 

0.369 

0.270 

0.503 

0.709 

1.173 

30 

0.232  . 

148/. 256 

0.569 

1.156 

32 

0.1J.6 


0.465 

0.691 

0.918 

TAM.L 


f  (cps) 
k 

=  0.5 
=  0.0140 

1.0 

0.C281 

8.0 

0.05'' 

4.0 

:  O.II25 

8.0 

0.2250 

18.0 

0.3375 

16.0 

0.4500 

aM° 

Pitching 

Amplitude , 

a  =  8° 

0 

0.127 

0.266 

0 . 468 

O.619 

3 

0.139 

0.286 

0.444 

6 

0.139 

0.296 

0.458 

9 

0.152 

0.310 

0.477 

12 

0.020 

0.059 

0.292 

-0.062 

-0.093 

0.248 

0.536 

15 

0.047 

0.034 

0.015 

0.052 

-0.221 

-0.098 

0.092 

18 

0.047 

0.170 

0.237 

0 . 190 

-0.258 

-O.606 

-O.302 

21 

0.066 

0.137 

0.294 

O.369 

-0.183 

-0.660 

-0 . 8o4 

24 

0.051 

0.141 

0.399 

0.430 

0.417 

-0.186 

-0 . 784 

27 

0.111 

-0.007 

0.330 

0.150 

0.595 

0.770 

1.028 

30 

0.085 

0.093 

0.508 

0-023 

0.458 

0.585 

1.085 

32 

0.028 

0.036 

-0.103 

0.145 

0.469 

0.802 

TABLE  II  ORIGINAL  TWO-DIMENSIONAL  AERODYNAMIC 
DAMPING  DATA,  HQp  ,  DOR  M  =  0.5 


f(cpa) 

K 

=  0.5 
=  0.0094 

1.0 

0.0187 

2.0 

0.0375 

4.0 

0.0750 

8.0 

0 . 150 

12.0 

0.225 

16.0 

0.300 

aM° 

Ditching 

Amplitude 

,  a  -4° 

0 

O.063 

0 . 176 

0.281 

0.392 

3 

O.O85 

0.183 

0.314 

0 . 36  6 

6 

O.O98 

0.209 

0.314 

0-399 

9 

O.O98 

0.222 

0.340 

0.444 

12 

0.059 

-0.052 

0.013 

0.092 

-0.150 

-0.013 

-0.477 

15 

0.007 

0.177 

0.320 

0.085 

0.288 

0.026 

-1.373/-. 889 

18 

O.O59 

0.601 

0.726 

.196/ -.281 

21 

1.431 

0.902 

1.712/. 758 

.431/. 706 

24 

O.863 

0.614 

0-595 

1.373/1.575 

27 

0.523  1.229/1.503 

1.203/1.144 

ClM0 

Pitching 

Amplitude 

,  a  =6° 

0 

O.O96 

0.177 

0.285 

0.383 

3 

0.087 

0.168 

0.279 

0 . 398 

6 

O.O99 

0.186 

0.282 

0.412 

9 

0.038 

0.166 

0.311 

0.430 

12 

0.070 

0.215 

0.017 

-0.073 

-0.171 

15 

-0.058 

0.157 

0.134 

0.322 

0.180 

-0.523 

-0.250 

18 

o.uoy 

0.122 

0.261 

0.360 

0.317 

-0 . 1.10 

0-099 

21 

0.601 

0.424 

-0.032 

0.613 

24 

0.221 

0 . 749 

0.651 

27 

0.273 

0.468 

0.511 

0.186 

30 

0.192 

0.325 

0.314 

O.58I 

32 

0.131 

0.224 

0.186 

O.58I 


ThbLK  IT  -  CONCLUDE!' 


f(cps) =  0.5 

k  =  0 . 0094 

1.0 

0.0187 

2.0 

0.0375 

4.0 

0.0750 

8.0 

0.150 

12.0 

0.225 

16.0 

0 . 300 

0  M° 

Pitching  Ampi t,ud 

3,  a  =  8° 

0 

0 . 0&  j 

0.193 

0.312 

0.407 

3 

0.085 

0.186 

0.314 

6 

0.093 

0.214 

0.30Q 

9 

0.131 

0.021 

0.227 

12  -0.013 

-0.087 

O.O95 

-0.154 

0.240 

-0.049 

0.291 

15  -0.025 

-0 . 103 

-0.134 

0.294 

•358/. 173 

0.322 

0.103 

18  -0.0u2 

-0.003 

0.172 

0.324 

0.387 

15 7/ -.229 

.002/. 060 

21  0.078 

0.047 

0.245 

0.358 

0.485 

O.I63 

0.216 

24 

0.221 

-0.526 

-0.371 

27  -0 . 054 

-0.002 

-0 . 198 

0.340 

0.449 

0.544 

30  -0.075 

-0.162 

-0.0o2 

0.201 

0.216 

O.567 

TABLE  III 

ORIGINAL  TV 0 -1) IMEN S I ONAL 
DAMPING  DATA,  Hq2,  FOR 

AERODYNAMIC 
M  =  0.4 

f(cps)  = 

0.5 

1.0 

2.0 

4.0 

8.0 

12.0 

16.0 

0.007 

0.014 

0.028 

0.0562 

0.1125 

0 . 1688 

0.225 

Pitching  Amp.,  itude 

,5=4° 

0 

0.072 

0.124 

0.222 

O.307 

3 

0.072 

0.144 

0.216 

0.314 

6 

0.085 

0.163 

0.242 

0.320 

9 

0.144  . 

150/. 137 

0.281 

.373/. 281 

12 

0.013 

0.026 

0.098 

0 . 183 

0.412 

0.484 

-.170/. Ill 

15 

0.059 

-0.C72 

O.I96 

0.039 

0.418 

0.543 

0.569 

16 

0.320 

0.575 

1.013 

0.771 

21 

0.373 

0.131 

0.765 

1.020 

24 

O.I83 

0.523 

1.020 

aM° 

Pitching  Amplitude,  3  =  6° 

0 

0.064 

0.122 

0.206 

0.288 

3 

0.061 

0.139 

0.212 

0.296 

6 

0.093 

0.186 

0.267 

O.369 

9 

0.090 

0.174 

O.183 

.I63/. 232 

12 

0.006 

0.020 

0.006 

0.134 

0.116 

O.163 

0.102 

15 

0.015 

0.023 

0.110 

0.078 

0.038 

0.067 

18 

0.055 

-0.023 

0.041 

0.151 

0 . 389 

0.177 

0.090 

21 

0.343 

0.378 

0-392 

0. 305 

2k 

0.035 

0.256 

Qm° 

Pitching  Amplitude 

,  a  =  8° 

0 

O.060 

0.136 

0.214 

0.320 

3 

0.065 

0.149 

0.239 

6 

0.098 

0.126 

0 . 147 

a 

0.II3 

0.14? 

0 . 14 1 
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Figure  1.  Two-Dimensional  Channel  for  UAC 
8-Foot  Octagonal  Wind  Tunnel. 

U2 


Figure  3.  Oscillograph  Traces  of  Pressures  on  Oscillating 
Two-Dimensional  NACA  001.2  Airfoil;  M  =  4,  f  =4cps. 
c  -  2  it,  k  =  057,  a  =  6°. 
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Figure  5-  Unsteady  Pressure  Coefficient  Magnitude  Versus 
Chord  Station;  a  =  6°,  M^.3. 
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Figure  9-  Effect  of  Kean  Incidence  Angle  on  Moment 

Hysteresis  Loops;  q  =  6°,  M  = .  2  ,f  =  4  cps  ,  k  =  .l  1 2. 
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REDUCED  FREQUENCY  PARAMETER,  k 


Figure  13- 


Var  i  at  i  on 
Parameter 
Incidence 


of  Two-Dimensional  Aerodynamic  Damping 
With  Reduced  Frequency  and  Mean 
Angle;  a- 6°,  M  =  .2. 


dimensional  aerodynamic  damping  parameter. 
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Figure  16. 


Representative  Cross  Plots  of  Smoothed  Two- 
Dimensional  Aerodynamic  Damping  Parameter. 
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Figure  17*  Two-Dimensional  Aerodynamic  Damping 
Surface  for  M  =  0.2. 
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Figure  19.  Two-Dimensional  Aerodynamic  Damping 
Surface  for  M  =  Q.U. 


Figure  20.  Comparison  of  Present  Data  and 
Reference  1  Data. 
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Figure  21.  Comparison  of  Present  Data  and 
Reference  1  Data. 
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TWO-DIMENSIONAL  AERODYNAMIC 
DAMPING  PARAMETER,  Ha2 


STALL  ANGLE  PARAMETER,  a. 


Comparison  or  present  Data  and 
Reference  lit  Data. 
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Figure  2h.  Variation  of  Aerodynamic  Damping  With 
Azimuth  for  an  S-61F  Rotor  Blade. 
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Figure  25-  Variation  of  Aerodynamic  Damping  With 
Azimuth  for  an  S-61F  Rotor  Blade. 
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Figure  2 6.  Effect  of  Incidence  Angle  Changes  on 
Rotor  Stability. 
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Ugo.re  27.  Reduced  P’requency -Incidence  Angle 
Trajectories  for-  Various  Loading 
Conditions  on  an  B-6l  Rotor  Blade. 
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STALL  FLUTTER  FLIGHT  CONDI TI OK  BOUNDARIES 


PURPOSES  AND  METHODS 


Upon  completion  of  the  stall  flutter  analysis  described  previously,  it  was 
recognized  that  stall  flutter  stability  analyses  for  a  more  comprehensive 
set  of  rotor  loading  conditions  were  needed  to  provide  some  definition  of 
those  flight  conditions  which  would  be  particularly  troublesome  from  the 
standpoint  of  stall  flutter. 

The  S-6l  rotor  was  chosen  as  a  medium  for  those  additional  studies,  and 
the  Normal  Mode  Transient  Analysis  was  used  to  calculate  required  input 
data  for  the  stability  analysis.  In  tor  r  utat  i  al  tip  speed  was  kept  at 
660  ft/sec.  Rotor  shaft  angle  war  1 al ] y  zero,  and  uniform  inflow 

was  assumed.  It  should  be  noted  that  elastic  blade  deflections  were  in¬ 
cluded  in  the  analyses.  Rotor  flapping  was  trimmed  within  2  degrees  for 
all  cases,  while  collective  pitch  was  changed  to  produce  variations  in 
rotor  lift  and  torque.  Advance  ratios  of  .[>17,  .1*73,  .422,  .381,  .317,  and 
.256  were  chosen,  corresponding  to  forward  speeds  of  210,  185,  165,  l*+9, 

124,  and  100  knots  respectively.  At  each  advance  ratio,  from  five  to  eight 
collective  pitch  settings  were  found  which  caused  the  rotor  blade  to  be 
torsionally  unstable  for  part  of  the  azimuth,  in  addition  to  one  or  two 
which  caused  it  to  he  completely  stable.  The  Normal  Mode  Transient  anal¬ 
ysis  provided  azimuthal  and  radial  variations  in  angle  of  the  attack  and 
Mach  number  for  input  to  the  stall  flutter  stability  analysis  described 
earlier.  These  data  included  elastic  blade  deformations  due  to  inertial 
and  quasi-steady  aerodynamic  forces,  including  stall  and  compressibility 
effects.  The  stability  analysis  program  provided  data  such  as  that  shown 
in  Figures  24  through  26.  The  Sikorsky  Aircraft  test  data  for  Mach  number 
M  =  0.2  were  used,  since  these  data  were  the  most  comprehensive  and  appli¬ 
cable.  With  reference  again  to  Figure  25  as  an  example,  the  total  number  of 
degrees  of  azimuth  for  which  negative  three-dimensional  aerodynamic  damping 
existed  was  determined  and  designated  <gSF  .  This  was  used  as  a  parameter 
for  measuring  the  intensity  of  stall  flutter.  With  the  present  analysis, 
this  would  seem  to  be  the  most  suitable  parameter  for  judging  the  severity 
of  stall  flutter. 

RESULTS 


The  total  unstable  stall  flutter  azimuth  range  ySF  was  plotted  against 
rotor  lift  coefficient  solidity  ratio  Ci_/<t  for  the  advance  ratios  .537, 
.473,  and  .422.  These  plots  are  shown  in  Figure  28.  Note  that  these 
graphs  define  the  parameter  1 as  a  single-valued  function  of  Cl/ct 
for  each  advance  ratio.  Unfortunately,  similar  relationships  were  not 
generally  obtained  for  the  lower  three  advance  ratios.  It  was  found  in¬ 
stead  that  a  range  of  values  of  were  possible  for  a.  given  value  of 

Cl/o’  .  For  this  reason,  plots  of  the  form  shown  in  Figure  28  could  not 
be  provided  for  the  advance  -atios  .381,  .317  and  .256.  However,  the  re¬ 
sults  shown  in  Figure  28  were  cross-plotted  as  shown  in  Figure  29,  to  pro¬ 
duce  boundaries  on  a  CL/cr  versus  /x  plane,  as  originally  intended. 


In  order  to  produce  flight  condition  boundaries  which  would  include  the 
lower  advance  ratios,  some  other  rotor  performance  parameter  was  needed, 
and  the  rotor  torque  coefficient-solidity  ratio  was  chosen  for  this  pur¬ 
pose.  Figure  30  was  therefore  prepared  from  the  previously  calculated 
data.  The  corresponding  cross  plot  is  shown  in  Figure  31. 

DISCUSSION 


One  of  the  basic  aspects  of  flutter  on  helicopter  blades  in  forward  flight 
is  its  dependence  on  quantities  which  change  with  azimuth  angle.  This  is 
particularly  true  of  stall  flutter,  which  tends  to  occur  on  the  retreating 
blade  side  of  the  rotor  disc,  since  angles  of  attack  beyond  stall  occur 
there.  Thus,  as  would  be  expected,  the  stall  flutter  analysis  predicts  a 
limited  azimuthal  range  of  instability.  It  would  be  expected  in  practice 
that  blade  torsional  vibrations  would  increase  while  in  the  unstable  re¬ 
gion  and  decrease  or  damp  out  while  in  the  stable  azimuthal  region.  There¬ 
fore, it  is  obvious  to  expect  that  the  torsional  amplitude  reached  and  the 
number  of  cycles  of  significant  torsional  amplitudes  would  grow  with  the 
azimuthal  range  of  instability.  This  range  of  instability  is  conveniently 
obtained  from  the  present  stall  flutter  analysis.  More  elaborate  param¬ 
eters  might  be  proposed  for  use  with  the  data  in  their  present  form,  such 
as  an  integral  of  work  done  on  the  blade  by  the  aerodynamic  torsional 
couples  as  the  blade  vibrates  continuously  at  some  nominal  amplitude  and 
traverses  the  region  or  regions  of  instability.  In  view  of  the  limited 
range  of  data  available,  the  selection  of  the  extent  of  the  total  azimuthal 
range  of  instability  as  a  parameter  is  appropriate. 

The  curves  on  Figure  29  define  values  of  rotor  lift  coefficient  for  which 
the  stability  analysis  predicts  various  amounts  of  unstable  azimuth  range. 
At  a  given  advance  ratio,  increasing  rotor  lift  from  the  stable  region 
first  causes  a  small  region  of  instability  over  a  range  of  azimuth  angles 
on  the  retreating  blade.  This  area  increases  in  extent  with  rotor  lift, 
until  the  blade  angles  of  attack  reach  high  enough  values  to  enter  another 
area  of  stability.  The  region  of  instability  then  splits  into  two  smaller 
regious,  which  diminish  in  size  and  grow  farther  apart  in  azimuth  as  rotor 
lift  increases.  These  observations  are  illustrated  by  Figures  2b  through 
27  and  in  the  corresponding  earlier  discussion.  Note  that  the  boundary 
segments  which  are  below  the  maximum  value  of  <#<SF  at  each  advance  ratio 
on  Figure  29  correspond  generally  to  a  continuous  region  of  instability, 
and  will  probably  cause  the  highest  torsional  amplitudes. 

It  is  of  particular  interest  to  compare  Figure  29  with  a  chart  of  a  similar 
type  appearing  as  Figure  12  in  Reference  2.  The  slender  region  «^SF  =  100° 
on  Figure  29  approximates  an  extrapolation  of  the  points  plotted  for  the 
full-scale  H-21  test  on  Figure  12  of  Reference  2. 

The  practical  significance  of  the  boundaries  on  Figure  29  will  be  discussed 
next.  Since  the  S-6l  rotor  torsional  natural  frequency  and  the  expected 
stall  flutter  frequency  are  about  8  cycles  per  revolution,  each  period  of 
torsional  vibration  at  its  normal  tip  speed  is  represented  by  about  1*5 
degrees  of  azimuth.  Some  judgement  is  required  in  order  to  stipulate  how 
many  cycles  of  instability  are  required  to  produce  troublesome  amounts  or 
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torsional  vibration.  Comparison  with  the  above  results  of  Reference  2 

1  1  _  ,  ,  1  ,  _  1  ,  _  ,  , _ _  ,  1  ^  J  a.  ,  ,  nl>J  1  «  W  ill  —  Q  0 

vouxu  tnuv  at  xeas  c  two  t,vtica  ui  uiamLnii^,  w,  ysr  —  , 

are  required. 

The  largest  values  of  yjr  obtained  at  the  lower  three  advance  ratios 
(  fx  =  .331,  .317 ,  and  .2^*0)  occurred  at  values  c*  Cl/o-  ,  which  were 
practically  coincident  with  the  test  values  for  the  11-21  rotor  on  Figure 
10,  Reference  2,  which  was  mentioned  curlier.  The  values  of  V'sr  were 
120°,  130°,  and  lko°,  respect  ively .  It  was  found,  however,  that  ysr 
could  vary  widely  with  little  or  no  change  in  Cl/<t  because  of  heavy  rotor 
stalling.  Thus,  the  trimmed  rotor  could  he  operated  at  advance  ratios  and 
lift  coefficients  close  to  those  given  above  with  very  little  stall  flutter 
instability.  The  occurrence  of  stall  flutter  was  found  to  be  more  the  re¬ 
sult  of  extensive  blade  stalling  than  the  result  of  attainment  of  a  partic¬ 
ular  value  of  rotor  lift. 

Since  the  occurrence  of  blade  stall  is  accompanied  by  \  rapid  rise  in  blade 
section  drag,  it  was  suspected  that  rotor  torque  would  provide  a  more 
definitive  indication  of  conditions  for  stall  flutter.  This  proved  to  be 
the  case,  and  Figures  30  and  31  show  the  resulting  boundaries.  It  will  be 
noted  that  the  lower  boundaries  occur  at  a  fairly  constant  torque  level  as 
advance  ratio  is  varied. 


The  calculation  of  a  few  rotor  loading  cases  for  the  untrimmed  rotor  and 
for  shaft  angles  other  than  zero  showed  that  the  stall  flutter  boundaries 
will  be  moderately  affected,  although  no  particular  trends  could  be  de¬ 
fined. 


The  radial  and  asimuthal  distribution  of  angle  of  attack  would,  of  course, 
be  somewhat  different  if  nonuni  form  inflow  were  used  in  place  of  the 
assumed  uniform  inflow.  Usually,  the  inclusion  of  nonuniform  inflow  tends 
to  reduce  angles  of  attack  on  the  retreating  blades;  hence,  the  boundaries 
presented  would  probably  be  shifted  upward  somewhat. 


&R  =  660  FT/SEC, 


Figure  29-  Stall  P’lutter  Boundaries  Related  to  Rotor 
Lift  Coefficient  and  Advance  Ratio;  S-6l 
Rotor,  12R  -660  r7/SLC,as  -0o,aiS  =  0o.  bls=0° 
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Figure  30.  Effect  of  Rotor  Torque  Coefficient  on  Stall 
Flutter  at  Various  Advance  Ratios;  S-6l 
Rotor,  fiR -660  FT/SEC,  as  =0°  ,a(S-  0°,  b,s=0° 
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Stall  Flutter  Boundaries  Related  to  Rotor 
Torque  Coefficient  and  Advance  Ratio;  S-6l 
Rotor,  ftR  - 660  FT/SEC,  as  =  0°.ols  =0°.t>i«^ 


CONCLUSIONS 


1.  The  measured  values  of  two-dimensional  aerodynamic  damping  parameter 

in  pitch,  Ba 2  ,  are  in  good  agreement  with  potential  flow  theory 

for  low  mean  incidence  angles,  aM  <  9°>  and  depart  from  the  theory 
for  aM  ^  12°. 

2.  At  high  incidence  angles,  a*  >  12°,  the  damping  parameter  in  pitch 

decreases  to  negative  or  unstable  levels;  it  then  increases  to  positive 
values  again  if  the  mean  incidence  angle  is  increased  enough.  This  is 
true  unless  the  reduced  frequency  k  is  smaller  than  about  .15. 

3.  The  Reference  1  data  and  the  present  data  are  in  generally  good  agree¬ 
ment  over  a  substantial  range  of  both  k  and  aM  .  Both  sets  of 
data  predict  essentially  the  same  two-dimensional  stability  boundary 
at  moderate  values  of  k  . 

4.  The  available  data  corresponding  to  M  =  0.2  are  adequate  for  the 
analysis  of  contemporary  rotor  systems .  More  extensive  ranges  of  re¬ 
duced  frequency  and  angle  of  ajbtack  data  are  needed  at  other  Mach 
numbers . 

5.  The  results  of  the  stability  analysis  for  the  S-6l  rotor  are  in  good 
agreement  with  those  of  Reference  2. 

6.  Use  of  the  Reference  1  data  yields  slightly  larger  blade  damping  for 
the  advancing  blade  but  produces  good  agreement  with  the  present  data 
for  the  retreating  blade,  even  when  instability  is  present. 

7.  The  rotor  lift  coefficient-3olidity  ratio  is  a  suitable  parameter  for 
definition  of  stall  flutter  boundaries  at  advance  ratios  of  .42  and 
above . 

8.  The  rotor  torque  coefficient-solidity  ratio  is  a  suitable  parameter 
for  definition  of  stall  flutter  boundaries  for  all  advance  ratios. 
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RECOMMENDATIONS 


1.  The  existing  data  descrihed  in  this  report  are  adequate  for  contem¬ 
porary  rotors  and  operating  conditions.  As  rotorcraft  forward  speeds 
increase,  it  is  evident  that  a  much  larger  range  of  data  will  be 
needed,  including  data  for  a  simulated  reverse- flow  region. 

2.  The  analysis  described  in  this  report  is  based  on  the  calculation  of 
torsional  work  for  a  hypothetical  cycle  of  torsional  vibration  which 
takes  place  at  a  number  of  discreet  azimuth  locations.  On  the  actual 
rotor,  the  azimuth  angle  and  therefore  the  mean  angle  of  attack  and 
relative  velocity  changes  continuously  during  the  torsional  vibration 
cycle.  In  future  stall  flutter  analysis  and  testing,  consideration 
should  be  given  to  the  time  history  of  the  large  variations  in  relative 
velocity  and  mean  angle  of  attack  that  are  actually  taking  place  during 
a  cycle  of  torsional  vibration. 
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appendix  I 

TWO  DIMENSIONAL  UNSTEADY  MOMENT  FORMULATION 


In  the  following  analysis,  it  is  assumed  that  the  two-dimensional  airfoil 
section  is  executing  a  single-degree-of-freedom  torsional  motion  about  a 
pivot  axis  at  a  .  This  motion  is  assumed  to  have  a  torsional  amplitude, 

5  ,  and  a  torsional  frequency,  u  =  2frf  .  (As  stated  earlier  in  the  main 

body  of  this  report,  complex  quantities  will  be  denoted  by  an  asterisk 
superscript,  and  amplitudes  of  sinusoidal  functions  will  be  denoted  as 
barred  quantities,  such  as  5  •) 

The  complex,  unsteady,  total  twisting  moment  M  *0T  experienced  by  the 
airfoil  about  its  pivot  axis  will  consist  of  the  sum  of  the  steady  mean 
moment  Mm  associated  with  the  mean  incidence  angle  <2m  and  the  unsteady 
moment  M*  associated  with  the  torsional  displacement  relative  to  the 

mean  incidence  angle  a* 


mtot  =  mm+  mu  (67) 

A  convenient  expression  for  the  unsteady  moment  due  to  a  single-degree-of- 
freedom  torsional  motion  about  a  pivot  axis  at  a  is  given  in  References 
8  and  13  as 

M*  =  7r^b4w2[Ma  -(La+Mh)(-i-  +  a)  +  Lh(-^  +a)2]  (68) 


In  Eq.  (68),  the  quantities  Lh  ,  La  ,  Mh  ,  Mq  8X6  unsteady  lift  and 
moment  functions  which  may  be  taken  from  any  convenient  theoretical  or 
empirical  source,  and  the  unsteady  torsional  displacement  has  been  written 
in  terms  of  the  torsional  amplitude  as 


*  s 

a  -  a  e 


(69) 


The  moment  coefficient  will  be  defined  by  the  general  equation 

M 


'M 


^PU2(2b)2 


and  after  a  bit  of  manipulation,  Eq.  (68)  becomes 

Cmu  =  ~^T~  [Mq  "(La+  Mh)(y  +  a)  +  Lh  (y  +  °) 


afiwt 


(70) 


(71) 


80 


This  formulation  will  be  returned  to  later.  Initially,  though,  the  anal¬ 
ysis  will  proceed  in  a  more  general  fashion,  beginning  with  the  coefficient 
form  of  Eq.  (67)* 


(72) 


and  if  it  is  assumed  that  the  unsteady  moment  coefficient  is  a  periodic 
function,  it  may  be  written  as 


=  V  * 

""u 


*  ~iuit 


(73) 


where  the  amplitude  £m.j  is  still  regarded  as  a  complex  quantity  to 
account  for  the  phase  shift  between  torsional  displacement  and  moment  re¬ 
sponse.  Thus,  Eq.  (72)  becomes 


-  r  .  7  *  .Jwt 
'■'U.  +  Cij  Q 
"TOT  mU 


=  c, 


"M 


(Cmur+ iC«m) 


,iu»t 


(7fc) 


and  in  terms  of  the  notation  of  Eq.  (7l)»  the  real  and  imaginary  parts 
are  given  by 

(75) 

(76) 


t  2 

cmur  =  ~ 2~ [mor~  ^LaR+ +aj  +  LhR^y  +aj  ja 
ff"u«  '  ~z~  [Mai  -(Lai+ Mhi)(y +a) +  Lf,x  (-^  +  a)2ja 


It  is  convenient  at  this  point  to  derive  the  expression  which  was  used  to 
provide  the  unsteady  theoretical  moment  loops  shown  in  Figure  10.  To  do 
this,  the  condition  of  25#  chord  pivot  will  be  imposed  on  Eqs.  (7^),  (75) 
and  (76),  whereupon  the  mean  moment  vanishes,  Cmm  ,  as  does  the  factor 
1/2  +0=0-  Then,  substitution  of  Eqs.  (75)  and  (76)  into  Eq.  (7I)  and 
expansion  of  the  exponential  factor  yields  the  equation 

Cmtot  :  yk2(MQR  +  i  Mai)  (  cos  cut  +  ism  tut)  a  (77) 


Si 


which  has  a  real  part  given  by 


Cmtotr  ?  yk2(Ma«  coscut  -  Mai  sm  cut )  a  (78) 

The  time  variable  in  Eq.  (78)  may  be  eliminated  by  manipulating  the  ex¬ 
pression  for  the  total  incidence  angle,  which  consists  of  the  sum  of  the 
mean  incidence  angle,  <3m  »  and  the  torsional  displacement  relative  to  the 
mean  incidence  angle,  a *  ,  or,  by  using  Eq.  (69), 


aTOT  :  aM  +  a*:  aM  +  aelu'  (79) 

The  real  part  of  Eq.  (79)  is  given  by 

aR  =  aM  +  a  cos  cut  (80) 

which  may  be  solved  for  cos  cut  as 

cos  cu  t  =  °ft_— --M-  (  ^ 

a 

Use  of  the  well-known  trigonometric  identity  relating  sin2  cut  and  cos2 cut 
yields 

sin  cut  =  ±4-^/a2-(aR-aM)2  (32) 


and  substitution  of  Eqs.  (8l)  and  (82)  into  Eq.  (78)  leads  to  the  result 


,mtot* 


f  k2[MaR(aR-aM)±MaIv/a^(aR-aM)2  ] 


(33) 


which  is  valid  for  an  airfoil  oscillating  in  pitch  about  its  23%  chord. 
For  use  elsewhere  in  this  report  in  comparing  theory  with  experiment,  the 
function  Ma  may  be  specialized  to  the  case  of  an  incompressible,  poten¬ 
tial  flow  past  the  airfoil  (Reference  8) 


Ma 


(81) 


whereupon  Eq.  (83)  reduces  to  the  form 


( a  R  -  a  M)  ± 


(85) 


In  this  equation,  the  negative  sign  is  associated  with  increasing  <1r 
and  the  positive  sign  with  decreasing  CtR 
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APPENDIX  II 

DERIVATION  OF  THREE-DIMENSIONAL  DAMPING  EQUATION 


The  equations  necessary  to  determine  the  three-dimensional  aerodynamic 
damping  will  now  be  derived.  Use  will  be  made  of  the  formulation  pre¬ 
sented  on  page  210  of  Reference  8,  with  the  notation  altered  to  conform 
with  present  usage.  As  in  Reference  8,  the  development  begins  with  the 
relations  for  the  virtual  work  due  to  a  torsional  displacement. 


8w 


Mu^Sa  (r)  -[  j  f Q<r)8aT 


(86) 


where 


Mu  2 


*/ 


(87) 


is  the  unsteady  moment  per  unit  span  and 

a  (r)  =  fa(r)  aT 


(88) 


is  the  spanwise  twist  amplitude  distribution,  5T  is  the  twist  amplitude 
at  some  reference  station  (in  the  present  case,  the  reference  station  is 
at  the  blade  tip),  and  fa  O’)  is  the  normalized  mode  shape  distribution. 
The  subscript  2  has  been  appended  to  the  unsteady  two-dimensional  moment 
to  distinguish  it  from  the  three-dimensional  quantity  which  will  be  in¬ 
troduced  presently.  The  generalized  force  per  unit  span  in  torsion  is 
obtained  by  dividing  the  virtual  work  by  the  virtual  displacement. 


Qq 


*/=  M.  /  Muz*)f 


SaT 


rr~  ro 


ttr) 


(89) 


Hence,  the  generalized  force  over  the  entire  span  (i.e.,  the  three-dimen¬ 
sional  unsteady  moment)  is  obtained  by  integrating  Eq.  (89)  over  the  span 
of  the  rotor  blade,  or 


Mu 


3*=X0V*  -7^fX?<a<r)a 


(90) 


Sit 


As  in  Eq.  (69),  the  two-  and  three-dimensional  moment  coefficients  are 
defined  by  the  formulas 


M 


U2 


=  ( r )  =  -^/3{r)u2(r )  [2b(r)]2  CMuJ(r) 


Muj  r  Jp  T  uT2  (2bT)2CMu* 


(91) 


where  the  subscript  T  denotes  values  taken  at  some  reference  station. 
After  Eqs.  (91)  are  substituted  into  Eq.  (90)  and  the  equation  is  divided 
through  by  the  reference  quantities,  the  result  is  given  by 


The  imaginary  part  of  each  side  of  the  equation  will  now  be  taken,  in 
amplitude  form,  and  use  will  also  be  made  of  the  dimensionless  spanwise 
variable 


1 


_  r~  rO 
'  rT-r0 


(93) 


whereupon  Eq.  (92)  becomes 


Cmu3i  = 


C*\j2Ifa(77)tf77 


(9M 


The  two-dimensional  moment  coefficient  amplitude  function  on  the  right-hand 
side  of  Eq.  (91*)  may  now  be  replaced  by  Eq.  (36);  and  after  using  Eqs.  (39) 
and  (88)  for  jx  and  01(77)  ,  respectively,  Eq.  (9U)  becomes 


CmU3I 


(t;  )  d  T7 


(95) 


As  In  Eq.  (35) ,  the  three-dimensional  aerodynamic  damping  in  pitch  vill  be 
defined  by  the  derivative  form, 


<— <a 


3 


da  t 


(96) 


After  the  operation  defined  by  Eq.  (96)  is  applied  to  Eq.  (95)  the  result 

is 


r'  (  P u2b2  \  ? 

~Q3  2  J0  L  uz  b2  J  aQ2  fa  m  dr!  (97) 


In  the  most  general  case,  Eq.  (97)  takes  into  account  the  possible  span- 
vise  variations  in  the  parameters  p  ,  U  ,  and  b  .  However,  in  the 
present  study,  both  the  density  and  the  semichord  vill  be  assumed  to  be 
constant  over  the  entire  active  span  of  the  rotor  blade,  and  only  the 
velocity  vill  be  considered  to  be  a  variable  quantity.  Thus,  after  de¬ 
fining  the  spanvise  velocity  ratio  by  the  formula 


V  (7))  z  \J{T))/\Jj 


(98) 


the  final  form  of  the  three-dimensional  damping  parameter  vill  be 


(99) 
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